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THE  GCARS  SYSTEM 

FORTRAN  IV  PROGRAMMERS  MANUAL 

PART  B  -  PROGRAMS  FOR  GENERATION  OF  ALTERNATIVES 

By 
A.  Keith  Turner 

INTRODUCTION 

The  GCARS  System  generates  a  series  of  alternative  highway  alignments 
by  analyzing  a  series  of  "route-independent"  highway  location  factors. 
Basic  information  on  each  factor  is  supplied  as  a  "value  surface" 
representing  the  potential  costs  or  benefits  associated  with  the  particular 
factor.  The  programs  concerned  with  the  development  and  testing  of  these 
value  surfaces  are  described  in  Part  A  of  this  manual  (  5  ). 

The  generation  of  alternatives  involves  the  following  four  steps  (  3  ) : 

(1)  Convert  value  matrices  to  utility  matrices  (value  surfaces  are 
always  described  in  a  matrix  form  within  the  computer). 

(2)  If  multiple  factors  are  being  analyzed,  combine  the  appropriate 
utility  matrices,  each  multiplied  by  a  selected  weighting  factor, 
to  produce  a  combined  utility  matrix. 

(3)  Develop  a  network  of  links  (a  utility  network)  joining  adjacent 
modes  on  the  utility  matrix  so  that  routes  can  be  selected 
across  the  area. 

(k)     Generate  a  series  of  alternatives  by  repeated  minimum  path 

analysis  to  the  network. 
This  portion  of  the  Programmers  Manual  describes  the  computer  routines 
developed  to  perform  these  steps.  The  routines  are: 


* 
Number  in  parenthesis  refer  to  references, 
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Program  NETSUM  -  Accepts  value  matrices,  converts  them  to  utility 

matrices,  and  combines  them  in  specified  ratios  to 
produce  a  combined  utility  matrix.  Constructs  a 
utility  network  link  table  suitable  for  the  minimum 
path  algorithm. 
Subroutine  MPATH  -  Performs  the  repeated  minimum  path  analysis  to 

generate  a  specified  number  of  alternatives. 
Subroutine  PRTOUT  -  Output  routine,  prepares  printed  output  for 

NETSUM  and  MPATH. 
Subroutine  PATHS  *-  Produces  a  map  of  the  generated  alternatives  on 

the  printer. 
The  following  two  routines,  described  in  Part  C  of  this  manual,  are 
used  in  conjunction  with  the  above  routines : 

Subroutine  STATS  -  (use  required)  -  called  by  NETSUM  during  the 

transformation  of  values  to  utilities,  and  also 
by  PRTCONU. 
Subroutine  PRTCONU  -  (use  optional)  -  may  be  called  by  NETSUM  to 

produce  a  printer- contour  map  of  the  combined 
utility  surface. 

METHOD  OF  OPERATION 

Figure  1  is  a  simplified  flow  chart  showing  the  transfer  of  control 
between  the  routines  during  the  analysis. 

Program  NETSUM  exerts  the  primary  control.  All  control  cards , 
containing  option  specifications  defined  by  the  engineer,  are  read  by 
this  program.  It  then  calls  on  the  associated  subroutines  as  required. 
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PRESENT  CAPACITIES,  OPTIONS,  AMD  RESTRICTIONS 

a)  Size  of  Value  Matrices 

The  present  maximum  size  of  value  matrices  is  50  x  50.  Neither 
dimension  can  exceed  these  limits  without  correcting  the  DIMENSION 
statements  in  the  various  programs.  Changes  in  shape,  for  instance  to 
100  x  25,  can  be  accomplished  without  increasing  the  overall  storage 
requirements  of  the  programs.  An  overall  increase  in  size  will  require 
a  complete  redimensioning  of  most  variables. 

Special  Note;  The  value  matrices  (Z,  ZZ)  are  dimensioned 
Z(NR0W,  NCOL). 
where  NROW  =  number  of  rows  (height  of  matrix) 

NCOL  =  number  of  columns  (width  of  matrix) 

b)  Number  of  Links  in  Network 

The  maximum  number  is  98OO.  This  corresponds  to  the  number  required 
by  a  50  x  50  matrix  with  four  links  out  of  every  central  nodes,  three 
out  of  edge  nodes,  and  two  out  of  corner  nodes. 

c)  Spacing  of  Network 

The  present  networks  are  assumed  to  be  composed  of  square  grids  -  the 
spacing  between  columns  equals  the  spacing  between  rows.  This  can  be 
modified  by  changing  values  of  XDIST  and/or  YDIST  in  NETSUM  (lines  1&  and 
19  page  B-30)  to  reflect  the  relative  spacings  of  columns  and  rows. 
Further  corrections  would  be  required  to  correct  the  maps  produced  by 
PRTCONU  and  PATHS. 

d)  Number  of  Factors  Included  in  One  Analysis 

The  present  maximum  number  is  16.  This  can  be  changed  merely  by 
redimensioning  VAL  in  Program  NETSUM. 
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e)  Input  Options 

The  default  option  is  card  decks.  However  an  optional  request  allows 
the  value  surfaces  and  their  leader  cards  to  be  read  from  standard  GCARS 
binary  tapes.  The  characteristics  of  such  tapes  are  described  in  the 
Users  Manual  (  k    ).  A  suitable  tape  can  be  developed  by  Program 
LOADG  (  5  )• 

f)  Output  Options 

Standard  output  includes  (l)  table  of  factors  and  weights,  (2)  table 
of  alternatives  showing  path  totals,  relative  path  values,  and  number  of 
links  per  path  (see  page  B-55)»  and  (3)  map  of  paths  of  alternatives. 

Debug  option  gives  additional  output,  but  tends  to  destroy  neatness 
of  some  standard  output  tables.  Additional  items  include  (l)  listings  of 
input  value  matrices,  (2)  listings  of  transformed  utility  matrices, 
(3)  listing  of  combined  utility  matrix,  (k)   link  summaries  for  all  paths. 

Punch  option  produces  punched  deck  of  combined  utility  matrix  (see 
page  B-57). 

Contours  option  produces  printer  contour  map  of  combined  utility 
matrix.  Various  contour  displays  are  possible,  for  details  read  the 
section  describing  Subroutine  PRTCONU  in  Part  C  of  this  manual  (  6  ). 

FUTURE  M3DIFICATIOM3 
The  following  are  eight  options  which  might  be  considered  in 
applying  GCARS  to  actual  location  problems.  Many  of  these  modifications 
or  options  require  only  a  small  amount  of  new  programming.  The 
following  items  might  be  changed: - 

(1)  Size  and/or  shape  of  value  matrices. 

(2)  Shape  of  grid  (at  present  a  "square  grid  is  assumed). 
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(3)  Central  percentage  of  links  to  be  removed  from  a  developed 
path  before  generation  of  next  alternative  (at  present  93% 
of  links  removed). 

(h)  Optional  complete  removal  of  certain  links  or  nodes  from  a 
network  to  completely  prevent  travel  over  certain  subareas. 

(5)  Retain  links  in  network  over  certain  control  sections  so  that 
many  choices  can  utilize  a  particular  bridge  or  tunnel 
crossing  site. 

(6)  Increase  the  maximum  number  of  choices  to  be  generated 
(present  maximum  =  7). 

(7)  Utilize  table  of  nodes  passed  through  by  each  alternative  to 
search  and  tabulate  values  of  various  factors  along  each  link 
of  each  choice. 

(8)  Revise  method  of  weighting  topographic  residuals  so  that  cuts 
are  no  longer  rated  equal  to  fills,  so  that  deep  cuts  or  fills 
are  weighted  more  expensive  per  foot  increase  after  an  initial 
shallow  height  or  depth. 

METHOD  OF  COMPUTATION 
a)  Program  NETSUM 

Figure  2  is  a  flowchart  for  Program  NETSUM.  After  reading  the 
control  cards,  NETSUM  reads  a  GCARS  header  card  and  the  value  matrix  (Z) 
from  either  tape  or  cards.  The  variable  ITYPE  on  the  header  card 
specifies  the  type  of  transformation  that  will  be  applied  to  convert  the 
values  to  utilities. 

If  ITYPE  =  1,  residual  values  are  to  be  minimized.  A  residual  of 
-2.0,  representing  a  fill  condition,  is  weighted  equal  to  a  residual 
of  +2.0,  representing  a  cut  of  equal  depth.  This  is  accomplished  by 
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FIGURE     2  .     FLOWCHART  FOR  PROGRAM  NETSUM. 


setting  all  values  of  the  Z  matrix  to  their  absolute  values;  Z  (l,j)  = 
ABS  (Z  (l,j)  ).  Then  Subroutine  STATS  is  called  to  determine  the 
maximum  and  minimum  values,  these  are  used  to  standardize  the  Z  matrix 
values  between  0.0  and  10.0  by: 

Z  (I, J)  =  ((Z(I,J)  -  MINZ)/(MAXZ-MINZ))*10.0. 
If  ITYPE  =  2,  the  values  are  in  terms  of  costs  which  only  need 
standardizing  to  a  range  of  0.0  to  10.0.  This  is  accomplished,  as 
above,  by  calling  Subroutine  STATS  to  determine  the  maximum  and  minimum 
values,  then  dividing  each  matrix  value  by  the  range  and  multiplying 
the  result  by  ten. 

If  ITYPE  =  3t  the  values  are  in  terms  of  benefits  and  must  therefore 
be  standardized  then  subtracted  from  ten  to  invert  the  "benefits"  to 
"costs". 

Z(K,J)  =  10.0  -  ((Z(l,j)-MENZ)/JttXZ-MENZ))*10.0 
The  resulting  utility  matrix  is  then  multiplied  by  its  appropriate 
weighting  factor  (VAL(n))  and  summed  with  all  preceding  utility-  weight 
products  to  produce  the  combined  utility  matrix  (ZZ). 

ZZ(I,J)  =  Z(I,J)*  VAL(N)  +  ZZ(I,J) 
In  the  above  equations : 

N  =  number  of  factor  in  sequence  of  factors 

I  =  Row  coordinate  of  element 

J  =  column  coordinate  of  element 

Z  =  utility  matrix  for  this  factor 
ZZ  =  combined  utility  matrix  (initially  set  =  0.0) 
When  completely  generated,  the  ZZ  matrix  is  transferred  to  the  Z 
matrix  locations  (z(l,j)  =  ZZ(l,j))  and  the  ZZ  locations  are  released  for 
other  uses. 
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NETSUM  allows  the  origin  and  destination  nodes  of  the  alternatives 
to  be  entered  by  row  and  column  coordinates.  These  must  be  converted  for 
the  minimum  path  algorithm  to  a  single  quantity  as  follows: 
origin  =  NHOME  =  (iX-l)*  NCOL  +  IY 

destination  =  NDEST  =  (JX-l)*   NCOL  +  JY 

Figure  3  defines  these  variables. 

The  link  table  is  next  constructed  from  the  combined  utility  matrix 
elements.  The  value  of  each  link  ( TLINK)  is  the  average  of  the  utilities 
of  its  two  end  points.  The  links  are  ordered  in  the  table  according  to 
their  origin  nodes;  the  first  origin  node  is  the  top  left-hand  corner 
element  (Z(l,l)),  the  last  origin  node  is  the  bottom  right  hand  corner 
(Z(NROW,NCOL)). 

Central  nodes  in  the  network  have  four  links  radiating  from  them. 
These  are  computed  and  stored  in  the  order  shown  in  Figure  k. 

Nodes  in  the  four  corners  of  the  network  can  have  only  two  connecting 
links,  those  on  the  sides  only  three.  However  on  large  matrices  the 
link  to  node  ratio  approaches  four.  As  adjacent  nodes  are  used  in 
sequence,  a  network  is  slowly  completed  that  includes  links  allowing 
travel  from  Nl  to  N2  and  N2  to  Nl,  in  other  words,  an  "undirected" 
network  (3  ). 

The  minimum  path  algorithm  requires  three  pieces  of  information  for 
each  link:--(l)  its  origin  node  (N),  (2)  its  destination  node  (j),  and 
(3)  its  utility  value  (TLINK).  The  original  version  of  the  algorithm  thus 
required  three  storage  locations  for  each  link  and  for  large  networks 
this  rapidly  filled  the  available  core  memoiy  on  the  CDC  6500  computer. 
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The  storage  requirements  were  therefore  reduced  by  "packing"  the 
origin  and  destination  information  into  a  single  number  called  NXJ 
according  to  the  simple  formula: 

NXJ  =  N*IPACK  +  J 
where     NXJ  =  combined  origin-destination  code. 
N  =  origin  node 
IPACK  =  a  very  large  constant  number 
J  =  destination  node 

The  selected  value  of  IPACK  was  UoOOOq  whose  decimal  equivalent  is 
I638U.  This  value  is  large  enough  to  prevent  computational  problems  in 
networks  joining  2500  nodes  (i.e.  50  x  50  matrix),  but  should  be  increased 
if  the  network  dimensions  are  increased.  The  HXJ  values  are  "unpacked" 
by  Subroutine  MPATH  to  obtain  the  appropriate  N  and  J  values,  through 
the  use  of  integer  arithmetic  as  follows: 
N  =  NXJ  /  IPACK 
J  =  NXJ  -  N*IPACK 

Figure  5  shows  a  sample  utility  matrix,  and  the  computed  utility 
network.  The  resulting  link  table  is  shown  in  Figure  6. 
b)  Subroutine  MPATH 

Subroutine  MPATH  utilizes  a  minimum  path  algorithm  developed  by  the 
British  Road  Research  Laboratory  (8  ).  This  RRL  algorithm  was  coded 
in  FORTRAN  by  Martin  (2  ).  MPATH  utilizes  a  modified  version  of  the 
Martin  coding;  the  modifications  are  summarized  in  a  following  section. 

Figure  9  is  a  flowchart  summarizing  the  operation  of  this  algorithm. 
Three  tables  are  used  in  a  computation,  a  TREE  table,  a  CUMULATIVE  table, 
and  the  prepared  LINK  table.  The  computations  required  to  determine  the 
minimum  path  from  node  12  to  node  1  on  the  sample  network  in  Figure  5 
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FIGURE     9.     FLOWCHART    FOR   SUBROUTINE    MPATH. 
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are  shown  in  Figure  7.  Figure  8  shows  the  stages  in  the  development  of 
the  partial  tree. 
i)  Initialization  Procedures: 

The  NCUM,  TCUM  and  NN  entries  in  the  CUMULATIVE  AND  TREE  tables  are 
set  to  zero.  The  TSUM  values  in  the  TREE  table  are  set  to  arbitrarily 
extreme  values  (99999-99).  The  TSUM  entry  for  the  origin  node  (NHOME) 
is  now  set  to  zero, 
ii)  Development  of  Partial  Tree: 

The  LINK  table  is  searched  to  find  an  entry  having  NHOME  as  an 
origin  node.  A  modification  was  made  to  the  original  Martin  FORTRAN 
coding  at  this  point.  Martin  arranged  for  the  search  of  the  LINK  table 
to  begin  at  the  nth  entry;  where  n  =  NHOME.   (in  the  example,  where 
NHOME  =  12,  the  search  would  begin  at  the  12th  entry).  This  will  save 
search  time  and  is  always  possible  since  the  link/node  ratio  is  greater 
than  one.  Note  however  in  the  example  network  (Figure  6)  the  link/node 
ratio  is  2.833  so  that  the  search  starting  at  entry  12  has  to  extend  to 
33  before  a  link  originating  at  12  is  encountered.  If  the  search  is 
begun  at  an  entry  some  multiple  of  NHOME,  according  to  the  link/node 
ratio,  the  length  of  search  is  greatly  reduced.  Accordingly  the  algorithm 
coding  has  been  modified  so  that  the  search  of  the  LINK  table  begins  at 
n  -  NHOME*LNKNDE,  where  ENKNDE  is  the  integer  value  of  the  link/node 
ratio  (in  the  example  LNKNDE  =  2).  The  search  now  begins  at  entry  2k, 
and  lasts  only  9  rather  than  21,  entries.  This  simple  modification 
speeded  up  the  computation  by  a  factor  of  two. 
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The  first  link  originating  at  NHOME  is  examined  and  the  destination 
node  J  is  determined.  If  the  TSUM  entry  for  J  is  99999-99,  J   has  not  been 
reached  yet  by  the  developing  tree.  The  link  is  therefore  entered  into 
the  CUMULATIVE  table,  TCUM  being  updated  to  reflect  the  cumulative  sum 
of  the  link  values  (TLINK)  from  NHOME  to  J.   If  the  TSUM  entry  for  J 
is  less  than  99999*99,  J  must  have  already  been  reached.  Since  the 
algorithm  always  finds  the  best  path  to  a  node  the  first  time  it  reaches 
it,  there  is  no  point  in  considering  other  links  leading  to  J.  Thus 
they  are  not  added  to  the  CUMULATIVE  table.  This  discarding  of  repetitive 
links  saves  computation  time  in  large  problems. 

The  NCUM  value  in  the  CUMULATIVE  table  is  set  to  reflect  the 
position  of  the  link  in  the  LINK  table.  It  thus  acts  as  a  pointer. 

The  next  entry  in  the  LINK  table  is  examined  to  see  if  it  also 
originates  from  NHOME.  If  it  does  its  destination  node  becomes  J,  and 
it  is  checked  and  added  to  the  CUMULATIVE  table  or  discarded  as  just 
described.  This  process  is  repeated  until  all  the  links  originating 
at  NHOME  are  found  (in  the  example  there  are  only  two). 

The  algorithm  then  picks  the  minimum  TCUM  value  in  the  CUMULATIVE 
table.  This  is  the  "cheapest"  path  from  NHOME  to  node  J  in  line  HCUM  of 
the  LINK  table.  Values  are  then  updated  in  the  TREE  table;  at  line  J  of 
this  table  NN  is  set  equal  to  the  line  number  in  the  LINK  table  and  TSUM 
is  set  equal  to  TCUM.  This  minimum  TCUM  entry  is  then  removed  from  the 
CUMULATIVE  table,  and  the  node  J  becomes  the  new  NHOME  (see  Figure  7). 

The  entire  above  process  is  repeated  with  this  new  NHOME,  further 
entries  are  added  to  the  CUMULATIVE  table,  and  the  minimum  TCUM 
is  selected  as  just  described.  The  process  continues  until  the  entry 
being  removed  has  a  J  which  corresponds  to  NDEST.  When  it  does  a  path  lias 
been  found  between  the  desired  origin  and  destination  so  that  the  process 
can  stop. 
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One  additional  check  is  made  during  the  operation  of  the  algorithm. 
As  each  minimum  TCUM  is  selected  from  the  CUMULATIVE  table  it  is  compared 
with  the  corresponding  TSUM  value  in  the  TREE  table,  to  see  if  it  is  still 
less.  Although  this  check  was  made  as  entries  were  being  placed  in  the 
CUMULATIVE  table,  it  is  possible  for  that  particular  node  to  have  been 
reached  prior  to  the  removal  of  this  particular  link  from  the  CUMULATIVE 
table.  If  this  has  happened  the  TSUM  value  will  be  less  than  the  TCUM 
value,  this  TCUM  value  is  therefore  eliminated  and  a  new  TCUM  searched 
for  this  occurs  twice  during  the  worked  example  (The  6th  and  7th  links 
being  removed  are  marked  with  a  *  to  show  this).  The  reader  is  invited 
to  work  through  the  sample  calculations  himself,  to  learn  the  operation 
of  the  algorithm, 
iii)  Development  of  Multiple  Alternatives: 

The  above  calculations  will  serve  to  locate  only  a  single  path  -  the 
"minimum"  or  "best"  path  -  between  any  origin- destination  node-pair. 
Further  modifications  are  required  to  develop  multiple,  distinct 
alternatives  between  an  origin  and  a  destination.  Ayad  considered  this 
problem  (  1  ).  Accordingly  the  algorithm  was  modified  to  reflect  his 
concepts.  The  process  is  quite  simple: 

1)  Find  the  minimum  path  between  the  origin  and  destination.  This 
becomes  a  "first  choice". 

2)  "Remove"  a  central  percentage  of  the  links  forming  this  path 
from  further  consideration.  Such  links  can  be  "removed"  by 
raising  their  TLINK  values  to  very  high  values,  making  them 
extremely  unattractive  during  further  studies. 

3)  Re-apply  the  minimum  path  algorithm  to  the  revised  network  to 
find  a  new  minimum  path.  This  will  be  not  quite  as  good  as  the 
first  minimum  path  located  earlier,  and  so  will  become  a  "second 
choice". 
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h)     Calculate  the  ratio  of  the  path  total  of  this  choice  to  the  path 
total  along  the  first  choice.   If  this  ratio  is  not  too  large, 

"remove"  the  central  litaks  of  this  choice  also  from  further 
consideration  and  return  to  step  (3). 
5)  Repeat  the  analysis  until  either  enough  alternatives  have  been 
generated,  or  until  the  ratio  of  the  currently  generated  path 
sum  to  the  first  choice  path  sum  exceeds  a  pre-determined  cut- 
off value. 
The  modifications  that  perform  this  process  are  coded  into  the  latter 
part  of  Subroutine  MPATH.  The  nodes  through  which  each  passes  are  stored 
in  a  matrix  called  NODE,  the  number  of  links  in  each  path  are  stored  in 
NUMLK.  These  values  are  utilized  by  Subroutine  PATHS  to  generate  a  map 
of  the  alternatives, 
iv)  Summary  of  Modifications  to  the  RRL  Algorithm: 

The  most  important  modifications  to  the  algorithm  are: 

(1)  changed  search  procedure  in  the  LINK  table,  which  speeds  up 
the  computations  by  a  factor  of  two. 

(2)  the  addition  of  methods  to  "remove"  links  from  the  network 
and  thus  generate  multiple  alternatives. 

(3)  the  storage  of  the  information  concerning  the  path  sums  (TOTAL), 
path  ratios  (RATIO),  number  of  links /path  (NUMLK),  and  the  nodes 
comprising  each  path  (NODE)  so  that  further  manipulations  are 
possible. 

c)  Subroutine  PATHS 

Subroutine  PATHS  produces  maps  of  the  generated  alternatives  on  the 
printer.  The  routine  begins  by  scaling  the  map;  given  an  overall  width 
for  the  final  map  in  inches  it  computes  the  maximum  number  of  print 
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positions  (characters)  it  can  assign  to  each  grid  across  the  page,  and 
how  many  rows  down  the  page.  TMs  is  accomplished  by  computing  the 
following: 

COL  =  WIDTH  *  10.0  +  2.0 
where     COL  is  an  integer  number  and  equals  the  number  of 

characters  allowable  since  there  are  ten  characters/inch. 

PS  =  (COL  -  00LS)/(C0LS  -l) 
where     COLS  =  cumber  of  columns  in  matrix. 

PS  now  is  the  number  of  characters  between  columns  (see  Figure  10). 

The  final  map  width  will  thus  be  COIMAX  =  FS*(C0LS-l)  +  COLS. 

It  will  never  be  greater  than  the  requested  size.  The  number 

of  characters  per  grid  is  thus  XSTOP  =  PS  +  1. 

The  number  of  rows  per  grid  (at  6  rows  of  type  per  inch)  is  thus 
YSTOP  =  ((PS+l)*6+5)/lO 

All  the  above  calculations  are  done  in  integer  arithmetic. 

The  actual  map  is  computed  and  printed  line  by  line,  to  save  storage 
requirements.  A  line  of  blanks  and  plus  signs  is  first  generated  (SYM)  to 
define  the  column  spacing  along  the  top  of  the  map.  Once  inside  the  map 
however  only  two  possibilities  exist 

l)  the  current  line  being  generated  passes  through  the  grid  points 


2)  the  current  line  is  between  grid  points. 

The  computations  also  reflect  the  fact  that  within  each  grid  there  is 
one  line  (the  first  or  top  line)  that  passes  through  the  grid  points, 
followed  by  a  sequence  of  one  or  more  lines,  each  of  which  will  be 
identical,  that  lie  between  grid  pd>ints.  This  occurs  because  the  paths 
travel  only  from  grid  point  to  grid  point  along  either  vertical  or 
horizontal  links  and  cannot  cut  across  grid  squares.  Figure  10  shows  the 
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defined  variables.  The  program  begins  wi th  a  line  of  blank  characters. 
The  program  next  calculates  if  the  line  is  a  top  line  through  grid  points 
or  an  intermediate  line.  If  it  is  a  top  line  it  is  only  necessary  to 
check  each  path  in  turn  to  see  if  it  passes  through  IZ  or  IZP.   If  it 
passes  through  either  one  or  the  other,  the  end  characters  should  be  set 
to  the  appropriate  number.  If  both  IZ  and  IZP  fall  on  the  same  path,  the 
path  must  pass  from  one  to  the  other,  and  all  intermediate  characters  are 
set  to  the  path  number. 

In  the  case  of  the  second  line  (the  first  intermediate  line)  all  four 
positions  IZ,  IZP,  IZDL  and  IZDR,  must  be  checked.   If  IZ  and  IZDL  both 
lie  on  the  same  path,  the  path  must  travel  up  the  left  edge  of  the  square. 
If  IZP  and  IZDR  both  lie  on  the  same  path,  the  path  must  travel  up  the 
right  edge  of  the  square.  When  either  of  these  situations  occurs,  the 
appropriate  end  character  is  set. 

This  analysis  is  repeated  for  each  square  along  a  line  and  for  each 
choice  starting  at  the  worst  choice  generated  (this  allows  better  choices 
to  superimpose  on  poorer  ones).  A  savings  in  computation  time  is  effected 
by  noting  that  the  third,  fourth,  etc.  lines  are  all  the  same  as  the 
second  line.  Once  it  is  generated,  it  can  be  printed  the  required  number 
of  times  without  recalculation. 

This  is  repeated  for  the  number  of  rows  in  the  matrix,  at  which  time 
a  map  of  all  the  alternatives  will  have  been  produced. 

PROGRAM  OPERATING  CHARACTERISTICS 
a)  Core  Storage  Required 

On  the  CDC  65OO  Program  NETSUM  and  its  associated  subroutines 
(including  PRTCONU  and  excluding  HYPSO)  can  be  compiled,  loaded,  and 
executed  in  ll6K  (octal)  storage  locations  when  dimensioned  to  handle 
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50  x  50  matrices  (maximum).  The  storage  requirements  will  be  reduced  if 
the  maximum  matrix  size  is  less.  For  instance  matrices  around  35  x  30 
can  be  handled  in  just  under  60K  storage  locations. 

b)  Running  Time 

This  is  somewhat  data  dependent,  and  thus  slightly  unpredictable.   I* 
is  linearly  dependent  on  the  number  of  choices  being  generated.  The 
average  time  required  to  generate  five  choices  for  a  combination  of 
several  factors  is  about  100  Sec  CP-time  and  15  Sec  PP-time.  This 
includes  compiling  time  as  well  as  the  production  of  a  contour  map  6f   the 
combined  utility  surface,  the  path  statistics,  and  a  map  of  the 
alternatives.  For  the  same  data,  by  way  of  comparison,  the  program 
compiled  and  two  choices  were  generated  with  similar  output  in  50  Sec 
CP-time  and  15  Sec  PB-time. 

c)  Control  Cards 

Figure  11  shows  a  typical  deck  set  up  for  the  generation  of 
alternatives.  The  control  card  parameters  are  defined  in  detail  on 
pages  B -28, B-2°»  under  the  description  of  NETSUM,  since  it  utilizes  all 
control  cards. 

The  first  two  control  cards  are  always  required.  The  contour  map 
control  cards  are  needed  only  when  the  contour  map  option  is  requested. 
The  header  cards  and  value  matrices  may  be  omitted  if  these  are  to  come 
from  tape.  Many  dfcta  sets  can  be  submitted  as  a  single  run. 

PROGRAM  LISTINGS  AND  DOCUMENTATION 
The  following  pages  contain  listings  of  the  various  routines  in 
order.  All  are  written  in  basic  FORTRAN  IV.   Each  routine  is  documented 
through  the  use  of  appropriate  comment  cards  at  the  beginning  and 
throughout  the  program. 
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Each  program  has  been  analyzed  by  Program  TIDY,  available  as  a 
library  program  at  Purdue,  which  resequences  all  statement  numbers, 
removes  extraneous  numbers,  and  generally  makes  the  program  neater.  For 
this  reason  the  statement  numbers  in  each  routine  increase  progressively 
in  increments  of  five,  and  all  FORMATS  (if  any)  are  collected  at  the  end. 
The  cards  are  each  given  a  sequence  number  which  appears  at  the  far  right. 
a)  The  Use  of  FORSTAT  Data 

Each  routine  was  next  analyzed  by  another  program,  FORSTAT.  This 
checks  the  syntax  of  source  (FORTRAN)  programs  and  produces  cross- 
reference  listings  of  all  variables  and  statement  numbers  found  within 
the  routine.  The  cross-reference  listings  utilize  a  series  of  "FORSTAT 
Numbers"  which  are  found  at  the  far  left  of  the  listing.  A  vertical  ink 
margin  line  has  been  added  to  separate  these  numbers  from  the  rest  of  the 
listing.  The  cross-reference  listings  are  found  at  the  end  of  each 
routine. 

All  symbols  defined  in  the  routine  are  listed  alphabetically  and  are 
identified  according  to  type,  definition  and  references  (where  referenced), 
The  following  symbols  are  utilized  in  the  cross-reference  listings :- 
Types         I  =  INTEGER 
R  =  REAL 

D  =  DOUBLE  PRECISION 
C  =  COMPLEX 
L  =  LOGICAL 
N  =  NAME  LIST 
ASF  =  ARITHMETIC  STATEMENT  FUNCTION 
Definition     C  =  COMMON 

D  =  DIMENSION 
E  =  EQUIVALENCE 
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All  other  terms  are  believed  self-explanatory.  These  cross- 
reference  listings  will  be  most  useful  should  any  reader  desire  to  make 
changes  to  these  routines. 
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b) 


LISTING    OF    PROGRAM    NETSUM 


PPOGRAM  NfT  SUM  (INPUT  .OUTPUT  ,  T  APE  1  »  HUNCH  t  TAPE5=I  NPUT  »  TAPE6=0UTPUT  ) 
PPOGRAM  NETSUM 


PURPOSE  

TO  PREPARE  VALUE  SURFACES  FOR  bF.NFRAlION  OF  ALTERNATIVES  .BY 
MINIMUM  PATH  ANALYSIS  (5UBHOUTINE  MPATH). 

NETSUM  TAKFS  ONE  OR  MORE  VALUE  MATRICES,  CONVERTS  THEM  TO 
UTILITY  MATRICES  BY  STANDARDIZATION  AND  INVERSION  (IF  NECESSARY). 
AN0  MULTTl-,  IES  THE,,,  hY  SUppL1tO  WEIGhTING  FACTriRS  To  PfioHUCE  THE 
COMRTNED  UTILITY  SUkFacE. 

1HE  COMBINED  UlILITY  SURFACE  IS  ANALYZED  TO  PRODUCE  A  LINK 
TABLE  . 

THE  GCARS  SYSTE~i 
PURDUE  UNIVERSITY 

A,  KEITH  TURNER. 
CIVIL  ENGINEERING, 
JUNE  ]96V. 


REOUIRFD  SUHROUTINtS 

-STATS 
-PRTOUT 
-MP4TH 
-PATHS 

OPTIONAL  SUBROUTINES — 

-PPTC0N4-  USED  TO  PRODUCE  PRIMER  CONTOUR  MAPS  Of  THE  COMBINED 

utility  surfaces.  (IF  mapping  option  requested.) 


INSCRIPTION    OF    CONTROL    CARDS    -- 

F IwST    CONTROL    CARD    -- 

IX-ICOLS  1-3) 
IY-(C0LS  4-6) 
JX-(COLS  7-9) 
JY-(COLS  in-12) 
NFAC-ICOLS  13-H) 

MCHCE-tCOLS  15-16)  - 

WIDTH-ICOLS  J7-2")  - 

MATRIX-(COI.S  ?l-2b)- 


IM-ICOI.S  26-35) 


IP-lCOLS  36-4(i) 


ITIX-ICOLS  U-*1') 


COLUMN  OF  DESIRED  ORIGIN  NODE 

ROW  OF  DESIRED  ORIGIN  NODE 

COLUMN  OF  DESlRtO  DESTINATION  NODE 

ROW  OF  DESIRED  DESTINATION  NOOE 

NUMBER  OF  FACTORS  TO  BE  INCLUDED  IN  THIS 

ANALYSIS 

NUMBFH  OF  ALTERNATIVFS  TO  BE  GENFRATFD 

(MAXIMUM  =  7,  DEFAULT  OPTION  =  5) 

WIDTH  OF  RATH  MAP  AND  CONTOUR  MAP  IN  INCHE 

(pUNCH  DECIMAL  PoInt»  maxImum  VALUp 

ALLOivtD  =  12.7,  DEFAULT  OPTION  =  1P.0) 

iht  word  *tape3*  indicates  the  value 
matrices  ahe  to  he  read  off  a  tape  rather 
than  from  cards,  thf  tape  should  contain 
A  standard  gcars  header  card,  followfo  by 
The  matrix  stored  row-wise, 
(default  option  - 

if  these  columns  are  blank,  thf  program 
will  expect  the  valuf  matrices  and  their 
header  carus  to  he  supplieo  as  data  cards. 

THE  wOkD  ^CONTOURS  *    CAUSES  SUBROUTINE 
PRTCON<,  TU  HE  CALLED  TO  PRODUCE  A  COMTOUR 
MAP  of  THE  FINAL  Com^'nED  UTILITY  SURFACE. 
(DEFAULT  OPTION-  NO  MAP.) 

THE  WORD  *PUNCH*  CAUSES  A  DECK  OF  CAPOS  TO 
BE  PUNCHED  CONTAINING  THF  FINAL  COMBINED 
UTILITY  SUHFACE. 
(DEFAULT  OPTION-  NO  CARDS.. 

THE  WORD  *U£HUG*  CAUSES  EXTRA  PRINTEn  OUT- 
PUT TO  BE  PRODUCED.   ADDITIONAL  OUTPUT 
INCLUDES  THE  INPUT  DATA  (VALUE)  MATRICES. 
THF  UTILITY  MATRICES,  AND  DETAILED  LINK 
SUmmARIES  oF  EACH  ALT£RNATIVF. 
(DEFAULT  OPTION-  STANDARD  QUTPllT  AN1  PATH 
MAP  ONLY. ) 


6 
B 

10 
12 
1* 
16 
1« 
2° 
22 
24 
26 
2fl 
32 
34 
36 
3B 
40 
4? 
14 
46 
»8 
50 
5? 
54 
56 
58 
6  0 
6? 
64 
66 
6H 
7p 
72 

7* 

7* 

78 

8* 

86 

BH 

9  0 

9? 

9* 

96 

9B 

10  0 

10? 

|Oi, 

106 

10B 

110 

HZ 

114 
116 
118 
120 
122 

i?6 

12B 
130 
132 
13* 
136 
13B 
1*0 
1*2 
I** 
1*6 
1*R 
,5" 
152 
15* 
156 
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C    SECOND  CONTROL  CARD  —  16* 


c 


166 


C  THIS  CARD  CONTAINS  *NFAC*  WEIGHTING  VALUES.   EACH  COMPONENT  16B 

C  UTILITY  MATRIX  WILL  BE  MULTIPItU  BY  THE  APPROPRIATE  WFIGHT  170 

C  HEFOHE  ADDITION  TO  FORM  THE  FINAL  COMBINED  UTILITY  MATRIX.  172 

C  THESE  WEIGHTING  FACTORS  MUST  THEREFORE  HE  IN  THE  CORRECT  ORDER,  174 

C  SO  THAT  THF  DESIRED  FACTOR  WILL  BE  MULTIPLIED  BY  THE  DESlPED  176 

C  WEIGHT.  178 

r  I80 

C  EACH  WEIGHTING  FACTOR  IS  ASSIGNED  A  FIVE  COLUMN  FIELD  <1-S»  182 

C  6-lnill-lSf  ETC.)  SO  THAT  UP  TO  16  FACTORS  CAN  BE  ENTERED  ON  A  184 

C  SINGLE  CARD.   IF  MORE  THAN  16  FACTORS  ARE  BEING  COMBINED  A  186 

C  SECon"  (OR  EV£n  THIRD)  CArO  SHOULD  HE  USED.   DECIMAL  POInTS  188 

C  SHOULU  HF  PUNCHED,  EACH  FACTOR  CAN  THEN  BE  PLACED  ANYWHERE  190 

C  WITHIN  THE  CORRECT  S  COLUMN  FIELO.  192 

C  (DFFAULT  OPTION-  NONE)  194 

C  l»6 

C    198 

C  200 

C    THE  FOLLOWING  CARDS  ARE  REOUIREO  ONLY  IF  A  CONTOUR  MAP  HAS  BEEN  202 

C    REQUESTED.-  204 

C  206 

C    - 208 

C  210 

C    MAP  PARAMETER  CARD —  212 

C  21* 

C  CON_(COLS  1-2)      -  THE  NUMBER  OF  CONTOUR  INTERVALS  (MAXIMUM  =  216 

C  19),   (  CON  =  DATA  RANGE  /  DESIRED  CONTOUR  218 

C  INTERVAL  ),  220 

C  T0UP-(COL  3>       -  IF  TOUR  =  o»  PROGRAM  USES  MAXIMUM  AND  222 

C  MINIMUM  FHOM  DATA  FOP  UPPER  AND  LOWER  CON-  224 

C  TOUR  LIMITS.  226 

c  IF  TOUR  =  1,  DESIRED  TOP  AND  BOTTOM  CON-  228 

C  TOUHS  REAU  In  FROM  SPECIAL  PROGRAM  CARO  230 

C  (SEE  BELow) .  232 

C  IF  TOUR  =  2,     VARIABLE  CONTOURS  RFAD  IN  234 

C  (LOW  TO  HIGH).   IN  THIS  CASt  SET  CON  EQUAL  236 

C  TO  THE  NUMBER  OF  VALUES  TO  BE  READ  IN.  238 

C  IF  TOUR  =  3,  CONTOUR  INTERVALS  BECOME  1/2  240 

C  STANDARD  DEVIATION  (MAX  =  *3  SIGMA,  MIN  =  242 

C  -1  SIGMA) .  2*4 

C  LINES-ICOL  4)       -  IF  LINES  =  0,  CONTOUR  LINES  ARE  PRINTED,  246 

C  IF  LINES  =  ,,  ALTERNATE  (EVEN  VALUFD)  348 

C  CONTOUR  BANDS  ARE  PRINTED,  250 

C  IF   .  INES  =  ?.%    A.L  BANDS  PRlNTED  (A,.  252 

C  PRINT  POSH  IONS)  ,  254 

C  256 

C    NEXT  CARD(S)  REQUIRED  ONLY  IF  TOUR  =  1  OR  2.  258 

C  260 

C  IF  TOUR=l,  THIS  CARD  SPECIFIES  MINIMUM  ANO  MAXIMUM  ELEVATION  262 

C  VALUES,  --  ZMIN(COLS  1-10)  AND  /MAX (COLS1 1-20 ) ,  264 

C  IF  TOUR=2,  THIS  CARO  SPECIFIES  A  FORMAT  (INCON)  TO  BE  USED  IN  266 

C  READING  IN  A  SERIES  OF  IRREGULARLY  SPACEO  CONTOUR  LINES.  THE  268 

C  VALUES  ARE  ARRANGED  Fr0M  L0WEsT  T0  HIGHEST  oN  SUBSEQUENT  CArDs.  270 

C  272 

C    2  74 

C    THE  FOLLOWING  DATA  CARDS  ARE  OMITTED  WHEN  *TAPE3*  INPUT  DATA  IS  276 

C    REQUESTED-  278 

C  -OTHERWISE-  280 

C  2B2 

C    THERE  SHOULD  PE  *NFAC*  DATA  SETS  EACH  CONTAINING  284 

C  1)   A  HEADER  CARO  WITH-  286 

C  -FACTOR-ICOLS  ]-4o)-  ANY  40  ALPHANUMERIC  CHARACTERS  TDENT-  288 

C  IFYING  FACTOR,  290 

C  -AREA-ICOLS  41-70)  -  ANY  30  ALPHANUMERIC  CHARACTERS  INDENT  292 

C  IFYING  AREA,  294 

C  -NCOL-ICOLS  71-73)   -  NUMBtR  OF  COLUMNS  IN  VALUE  MATRIX,  296 

C  -NROW-ICOLS  74-76)   -  NUMbtR  OF  ROWS  IN  VALUE  MATRIX.  298 

C  -ITYPE-ICOLS  77-78)-  IF  ITYPF  =  1,  VALUE  SURFACE  IS  TOP0-  300 

C  GRAPHIC  RESIDUALS.  302 

C  IF  ITYPF  =  2,  VAi_UE  SURFACE  MEASURED  304 

C  IN  TERMS  OF  COSTS.  306 

C  IF  ITYPE  =  3.  VALUE  SURFACE  MEASURED  308 

C  IN  TEHMS  OF  BENEFITS.  310 

C  2)   MATRIX  CAPOS  -  CONTAINING  VALUE  MATRIX  In  ROWS  IN  STANDARD  312 

C  FORMAT  (8F10.3).   EACH  NEW  ROW  STARTS  ON  A  FPESH  CARn,  AND  314 

C  IS  STORED  ON  NCOL/8  TAROS,  316 

C  318 

r 32„ 


B-30 


10 

c 
c 
c 


15 


c 

?0 


g5 
30 


C 

*0 


45 
50 


DIMENSION  10TALI7),  HAT10(7).  NUMLM7) 

DIMENSION  7(50,50).  ZZI50.50),  val(1*i 

nl„tN«-TON  F*CToH(4)i  AREA!  3) 

DIMENSION  IC0DE(4) 

DIMENSION  TITl  (8) 

OI"i£NSTON  NXJ(98(J0),  TLINM9R00),  UUMMYU614) 

DATA  (TITL  (I  >  »l  =  lib>/10MCONT<iUH  MA.10HP  OF  THE  C  ,  1  OHONSTRUC TFO , 1  OH 
1  UTILITY  S«  10HUHFACE     / 

DATA  (TCOOfc (I ) , 1=1,4) /ShPUnCh, lOMCPNTOUPS   ,  5H0E8UG, 5HTAPE3/ 

OATA  IPACK/40000H/ 

COMMON  NX J, TL i NK,/, 7  7, ARE  A, FACTOH»NCOL,NROw,KOUNT,NUMLK. K.N, J, PVAL 
i  ,  I'HJT  .VALUE  ,OUMMY»NMf)ME  ,NDE ST ,NCHCE  ,  IM,  IP,  IF  I  X  ,  NF  AC 

COMMON  /MAPS/  CON, TOllH, LINES. *  I D Th , I HYPSO » I NCON ( 8 )  ,URANG(2S>  ,?MIn, 
,  ZMAX 

INlEGER    CON,  tour 

INTEGER    SYSIN.SDF3 

DATA    SYSIN,SDF"3/5»3/ 

PF«L  MAXZ»M]N7 

nVAL (T, J) = < 1-1 )*nC()L*J 

XOIST=l »0 

TDTST»1 .0 


INPUT  KOUT1NF 


REAO  (SYSIN,140)  1A.IY. JX, JY.NFAC.NCHCE.WlDTH.MATRIX.lM.IP.IFIX 

IF  (FOF.SYS1N)  13b, 10 

READ  (SYSlN.145)   lVAL (I) ,I=1,NFAC) 

IF  (  IM.NE.ICODE  (?)  )  (SO  TO  25 


READ  CONTOUR  map  CONTROLS  ONLY  IF  CONTOUR  MAP  REQUESTED 


HEAD  (^YSIN.ISO)  CON, TOUP, LINES 
IHYPSO=0 

if  (took. GT,2. OR. Tour. le.O)  go  To  25 

GO  TO  (  15,20)  ,  TOUR 

kfcaO  (SY51N,155)  iMlNtZMAX 

GO  TO  ?5 

►<EA0  (SYS1N.160)  InCoN 

REAO  (SYS1N.INC0N)   (l)RANR(I  )  ,1  =  1  ,CUN) 
DO  30  l=j, 50 
DO  30  J=l,50 
ZZ(I,J) =0.0 
DO  10"  N=l,NFAC 
VA(  OF  =V«L  (FJ) 


RE  10  VAI  HE  MATRICES  FROM  CARDS  OR  TAPE  ACCORDING  TO  CONTROL  REQUEST 


IF  (MATRIX. EQ.ICOOE (4) )  GO  TO  40 

REAO  (SYSIN.165)  FACTOR, AREA. NCOL.NROW.ITYPE 

DO  3b  1=1 ,NROw 

REAO  (SYSIM.170)  (Z ( I, J) , J=1,nC0L) 

GO  TO  5(1 

REAU  (SOF3)  FACTOR.ARF A.NCOLtNHOW, ITYPE 

00  4  5  1=1 ,NROW 

READ  (SUFJ)   (7(1, J) ,J=] ,NC0L) 

IF      (N.Flj.  1)     CALL    PRTOUT     (1) 

CALL    PHTOIJT     (?) 


wRITF  DERUG  PRINTOUTS  ONLY  WHEN  REQUESTED 


IF  (IF  IX.EO.ICUDF ( 3) )  CALL  PRTOUT  (9) 
GO  TO  <5=.,7U»80)t  ITYPE 


CO'VFRT  VALUES  TO  UTILITIES  SO  AS  TO  MINIMIZE  RESIDUALS 


00  5  0  1=1 ,NROW 

00  fcO  J=l ,NCOL 

/ ( 1 , J)=AHS(7 ( t ,J) ) 

CALL  STAIS  (NM)w,NCOL.50,50.Z.*HAR,XP,ZSIG,MAXZ,MINZ) 

00  Fib  1  =  1  ,NROw 

DO  5b  J=l iNCOL 

/(I,J)=((7II . J)-MTNZ)/(MAXZ-MINZ) 1M0.0 

IF  <  IF  1X.F'  .KODF  (3)  >  CALL  PRTOUT  (10) 

GO  TO  MO 


32? 
324 

326 
328 
330 
332 
334 
336 
338 
3*0 
342 
344 
346 
348 
350 
352 
354 
356 
358 
360 
362 
364 
366 
368 
370 
372 
374 
376 
378 
380 
382 

384 
386 
368 
390 
392 
3«4 
396 
398 
400 
402 
404 
406 
408 
410 
412 
414 
416 
418 
420 
422 
424 
426 
428 
430 
432 
434 
436 
438 
440 
442 
444 
446 
448 
450 
452 
454 
456 
458 
460 
462 
464 
466 
468 
470 
47? 
474 
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CONVFhT  VuLUtS  TO  UlILIIItS  SO  AS  10  MINIMIZE  COSTS 


70 


7b 


r. 
c 
c 
c 

HO 


Hb 

c 


CALL  STATS  (NHOw , NCOL i50 »bO ,Z .ZbAR.XP , ZSIG tMAXZ ,MTNZ) 

00  7b  1=1 ,NROw 

00  7b  J=1»NC0L 

7<T,J)=< (Z(IfJ)-MINZ)/(MAXZ-MIN2) )*10.0 

IF  (IFIX.EB.ICUDE (3) )  CALL  PRTOUT  (10) 

GO  TO  9U 


CONVERT  VALUES  TO  UIILITIES  SO  AS  TO  MAXIMIZE  BENEFITS 


95 

100 


CALL  STATS  (NR0W,NC0L.50,50,Z,ZHAR,XP,ZSI6,MAXZ,MINZ) 

00  8b  1  =  1  ,NROW 

DO  85  J=1»NC0l 

Z( I, J) =10.0-( (Z( I,J)-MINZ)/ (MAXZ-MINZ) )»10.0 

IF  (IF  IX.EO.ICOOEO)  )  CALL  PRTOUT  (10) 


SUM  THE  COMPONENT  UIILITY  SURF ACE*wt IGHTS  TO  OBTAIN  FINAL  SUFACE 


DO  9b  1=1, NROw 

00  9b  J=l ,NCOL 

ZZ(T.J)=/(I,J)*VAL(N)*ZZ(I,J) 
CONTINUE 

NHOME*( IY-1) »NCOL*IX 
NUFST=(JY-1 ) »MC0L*JX 


END  OF  INPUT  ROUT  INF 


THIS  AREA  INCLUDES  CALLS  TO  PRTOUT  TO  LIST  FINAL  COMBINED  UTILITY 
SURFACE,  TO  PUNCH  OUT  MATRIX  IF  RtQUtSTtl) , AND  TO  MAKE  A  MAP  OF  THESU 


105 

IK 

c 
c 


IF  (IFIX.EO.ICODEO)  )  CALL  PRTOUT  (3) 

IF  (IP.EO.ICODE(l)  )  CALL  PRTOUT  (I*) 

IF  (T.M.NE.ICO0E  (?)  )  GO  TO  10*1 

CALL  PRTC0N4  (ZZ,UTL«NUOW,NCOL»50,50> 

00  110  1=1, NROW 

00  11(1  J=1,MC0L 

7<T,J)=ZZ(I,J) 


GENERATE  THE  LINK  TABLE 


NUM=1 

00  130  1=1, NROW 

IUPal-l 

IDNol+1 
00  13i  J=1,NC0L 

JLFT=J-1 

JRT=J*1 

IF  (JLFT.LT.l)  GO  TO  115 

NXJ(NUM) =  NVAL (I, J) »IPACK*NVAL (1,JLFT) 

TL I NK ( NUM >  =  <  Z  <  I .  JLF  T  >  ♦  Z  <  I ,  J )  >  « XO I  ST  /2  .  0 

Nt.JM  =  NUM+  ) 

IF  (JRT.GT.NCOL)  GO  TO  120 

NXJ(NUM)=NVAL(I. J) »IPACK*NVAL(1, JRT) 

TLINK(NUM)=(Z(1,JRT) ♦Z(I,J))«X01ST/?.0 

NUMsNUM* 1 

IF  (1UP.LT.1)  GO  TO  125 

NXJ(NUM)=NVAL(I, J)*IPACK»NVAL(1UP,J) 

TLTMK (NUM) =(Z(IUPiJ) *Z (I »J) )«XDI ST/2.0 

NUMsNUM* 1 

IF  ( ION. GT. NROW)  GO  TO  130 

NXJ(NUM)=NVAL< I,J)»IPACK*NVAL(1UN.J) 

TLINMNUM)  =  (Z<I[)N,J)*Z(I,J>  )«X01  ST/2.0 

NUM=NUM»1 

CONTINUE 

NSUMsNUM"! 


C    CALL  MINIMUM  PATH  ROUTINE   TO  GENERATE  ALTERNATIVES 


12" 


12b 


CALL  MPATH  (NSUM. WIDTH) 
GO  TO  5 


13s 


476 
478 
480 
482 
484 
486 
488 
490 
492 
494 
496 
498 
500 
502 
504 
506 
508 
510 
512 
514 
516 
518 
S2o 
522 
524 
526 
528 
53o 
532 
534 
536 
53H 
540 
542 
544 
546 
548 
550 
552 
554 
556 
558 
560 
562 
564 
566 
568 
570 
572 
574 
576 
578 
580 
582 
584 
586 
588 
590 
592 
594 
596 
598 
600 
60? 
604 
606 
60R 
610 
612 
614 
616 
61R 
620 
622 
624 
626 
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1X0. 

r 

14(1 

FORMAT 

<4I3.2I2tF4.0,A5,,A10»2A5> 

111. 

]4S 

FOHMAT 

(16F5.0) 

112. 

15" 

FORMAT 

(I2»2I1  ) 

113. 

155 

FOhMAT 

(2Fl0.ni 

114. 

160 

FORMAT 

(BAIO) 

115. 

165 

FORMAT 

(4A10»3A10»2I3»I2) 

116. 

l"Tn 

FORMAT 

(8F10.3) 

117. 

END 

626 
630 
632 
634 
636 
638 
640 
642 
644 
646- 


c)  FORSTAT    CROSS-REFERENCE    LISTINGS    FOR    PROGRAM    NETSUM 

PRCHiPnM  NFTSUM 

REFERENCES 


29. 


EM    NUMBER 

DEFINITION 

5 

20. 

10B. 

10 

22. 

21. 

15 

28. 

27. 

20 

30. 

27. 

?5 

32. 

23. 

26. 

30 

3'.. 

32. 

33. 

15 

40, 

39. 

40 

42. 

37. 

*5 

44. 

43. 

50 

45. 

41. 

5b 

49. 

46. 

6C 

51  . 

49. 

50. 

6K 

55. 

53. 

54. 

70 

58, 

46. 

75 

61  . 

59. 

60. 

B(J 

64. 

46. 

«fc 

67. 

65. 

66. 

4(1 

69, 

57. 

63. 

45 

71. 

69. 

70. 

10(1 

72. 

35. 

1    5 

7". 

It. 

IK 

H)  . 

79. 

80. 

lib 

93. 

89. 

12( 

97. 

93. 

125 

101, 

97. 

13(1 

105. 

83. 

86. 

135 

109. 

21. 

ur 

110. 

FORMAT 

20. 

U5 

111. 

FOHMAT 

22. 

150 

112. 

FORMAT 

24. 

155 

113. 

FORMAT 

28. 

160 

114. 

FOHMAT 

30. 

If  s 

115. 

FORMAT 

3b. 

170 

116. 

FOHMAT 

40. 

101, 


14    Sr/;TFV> 

NT    NUMBERS 

DEFINED 

3    M  (\  >  J  y  1 1 

■     NUMUth    OF 

RFFEKtNfES 

4  6    T  m  A  | 

NllMRfcP    OF 

RfcFKRKNCES 

1  .4    A,. r ••  Ai- 

NUMREH    Of 

RF.Ft  HI-  NirF.S 
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ai 
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d)  LISTING    OF   SUBROUTINE    MPATH 

1,  SUBROUTINE  MPATH  (M.INK, WIDTH) 

SUBROUTINE  MfATH 


PURPOSE  -- 


THE  SCARS  SYSTEM 
PURDUE  UNIVERSITY 


A.  KEITH  TURNER, 
CIVIL  ENGINEERING, 
JUNE  1969. 


MIN 
MIN 
-MIN 
MIN 
MIN 
MIN 

to  generate  a  series  of  alternatives  by  repeated  application  ofmin 

minimum  path  analysis  to  a  utility  nttwork  uevelopeo  by  program     min 

nftsum.  min 

this  subroutine  usts  a  modified  version  of  the  british  road    min 

research  laboratory  (rrl)  algorithm.  the  basis  of  this  fortran  coo-min 

ing  was  developed  by  -  min 

b.  v.  martin,  1963,  min 

minimum  path  algorithms  for  transportation  planningmin 

civil  engineering  systems  laboratory  reasearch     min 

report  r63-52,  min 

massachusetts  institute  of  technology,  min 

Cambridge,  mass.  min 

MIN 
................... MIN 

MIN 

MIN 
MIN 
MIN 

MIN 

MIN 

MlN 
MIN 

MIN 

MIN 
MIN 
MIN 
MIN 
MIN 
MIN 
MIN 
MIN 

..£ — MlN 
MIN 
MlN 
MIN 

TSUM  MIN 
MIN 
MIN 

—  .  —  MlN 
MIN 
MIN 
MIN 
MIN 
M.IN 
MIN 
MIN 
MIN 

.....MlN 

)      MlN 

(jOO.MIN 

mIn 

MIN 

MIN 

N, J.RMIN 

ifix,mIn 

MIN 
MIN 
MIN 
MIN 
MIN 
MIN 
MIN 
MIN 
MlN 
MIM 
MlN 
MIN 
MlN 
MIN 
MIN 
MIN 
MIN 


CALL  MPATH  (NLINK, WIDTH) 
WHERE  -  NLINK  =  NUMBER  OF  LINKS  IN  NETWORK 

_  WIDTH  =  WIDTH  OF  MAP  OF  PATHS  IN  INCHES  (TO  BE  PRO 
DUCED  BY  SUBROUTINE  PATHS) 


DESCRIPTION  OF  CONTROL  CARDS  — 

NONE  REQUIRED-  ALL  NECESSARY  VALUES  SUPPLIED  BY  PROGRAM  NE 

VI»  COMMON  STATEMENT 


NEOUIRED  SUBROUTINES- 
-P4THS 
-PRTOUT 

OPTIONAL  SUBROUTINES- 

-NONE 


DIMENSION  AREA (3),  FACT0W<4),  NXJ(980U),  TLINKI9B00),  NNI^sOO 
DIMENSION  TSUMIgbOO),  TCUMUOOO),  NCUM(IOOO),  ITABl^OO).  NODF 
17) 
OIMFNSION  TOTAL(7),  RATlOi?),  NUMLK(7),  ICODE,3) 
DATA  IPACK/40000B/, (ICODE (I) . 1=1 ,3 ) /SHPUNCH.5HMAP   »SHOEBUG/ 
COKMON  NXJ,TLINK,NN,TSUM,AHEA,FACTOR,NCOL,NROW,KOUNT,NUMLK,K, 

lVAl,InDT,VALUF,ToTAL,RATTo,ITAB,NODE,NHOME»NOEST,NCHCE,lMtIP, 
?NFAC 
DATA  (TCOOE(I) , 1=1,3) /5HPUNCH,5HMAP   ,5HDEBUG/ 


lNlTI«I.I7t  ALL  ARRAYS 


KOUNT=n 

N70NE=K'COL«NROW 

LNKNDF=NLINK/NZONE 

IF  (nChCE.lT.I.OR.NChCE.GT.7)  NCHCt=5 

CM I  PRTOUT  (111 

DO  10  I«1,NZ0NE 

TSUM(I) =99999.99 

MN (NHOmE) =0 

rSUM(NHOME) =0 

'MCMaO 

MMbNHOME 

N'MMbNH()ME*LNKNDE 


2 

4 

6 

8 

10 

12 

14 

16 

18 

20 

22 

24 

2ft 

28 

30 

32 

34 

36 

3H 

40 

42 

44 

46 

4H 

SO 

52 

54 

56 

58 

60 

62 

64 

66 

68 

70 

72 

74 

76 

78 

80 

82 

84 

86 

88 

90 

92 

94 

96 

98 

100 

102 

104 

106 

108 

110 

11? 

114 

116 

118 

120 

122 

12* 

126 

128 

130 

132 

134 

136 

138 

140 

142 

14* 

1*6 

1*8 

ISO 

152 

15* 

156 
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60 


SFAHCH    LINK    TABLE    FOR    LINKS    FROM    NtnUME  , 
ST»KT     SFAKCH    AT    FNTRY    NMM. 

IF     <  (NXJ(NMM)/IPACK)-NM)     25,20,20 
NHM=NMM- (LNKNTF*i  ) 
on    TO    15 

DO  4(1  T=NMM,M  INK 
NrNXJII  )  /IPACK 
IF  (N-NM)  40>30.45 
J=NXJ(I )-N«IPACK 
K  =  J 

IF  TSUM  OF  ENO  NODE  EQUALS  99999.99. 
ENTER  LINK  IN  CUMULATIVE  TABLE 

IF  (TSUM,K)-99999.99)  40, 3b, 40 
NC'M  =  NCM*  1 

TCUM(NCM) =TSUM(NM) *TLINK(I) 
NCUM (NCM) =T 
CONTINUE 

FIND  MINIMUM  TCUM  IN  CUMULATIVE  TABLE 

TMIN=99999.99 

DO  55  K=l ,NCM 

IF  (TMIN-TCUM(K) )  55.55,50 

TmTn=TCUm(K) 

l.eNCUM(K) 

MsK 
COMTlNUt 
N=NXJ(L)/IP'CK 
K=NXJ(L)-N»IPACK 
IF  (TSUM(K)-TMIN)  b0,60,65 
1  =  1 
(50  To  70 


C 

C    

C    FNTEP  THIS  MINIMUM  LINK  IN  TREE  TABLE 
C         . 

h5      TSUM(K)=TMIN 
NN(K)=l 

1=0 
c         

C    IF  END  MODE  EOUALS  DESTINATION,  PATH  FOUND.  GO  TO  90  . 
C    

IF  (N0FST-K)  70.9U.70 
C    . 

C    UPOATF  CUMULATIVE  TABLE  ENTRIES. 

C    N.AKE  LAST  tND  NODE  NHOmE  AND  RETURN  TO  15  . 

C    

70     00  75  NM=M,NCM 

TCUM (NM) =TCUM (NM+1 ) 
nCUm(Nm) =mCUm(Nm*1) 
IF  (NM+l-NCM)  /5.H0.75 
75     CONTINUE 
HO      NC"=NCM-1 

IF  (I)  45.8S, 45 
«5     NMMrK«LNKNDE 
NMrK 
OO  10  15 


DESTINATION  REACHED.  TRACE  PATH  FROM  TREE  TABLE.  STORE 
ICS. 


KOUNT=KOUNT*l 

IOUT=0 

IF  (If IX.EQ.ICODEU))  CALL  PRTOUT  (5) 

IF  (IFTX.FQ.ICODE (3) )  CALL  PRTOUT  (6) 

IOUT=IOUT*i 

NODE (lOUT.KOUNT)sK 

IF  (K-NHOME)   100.105,100 

L=mn(K1 

1  T^H  (  IOUT  )  =1 

K=''lXJ  (L)/IP»CK 

GO  TO  us 


MIN 
MlN 
MIN 
MIN 
MIN 
MIN 
MlN 
MIN 
MIN 
MIN 
MIN 
MIN 
MlN 
MIN 
MIN 
MIN 

mIn 

MIN 

MIN 
MIN 
MIN 
MIN 
MIN 
MlN 
MIN 
MIN 
MIN 
MIN 
MlN 
MIN 
MIN 
MIN 
MIN 
MlN 
MIN 
MIN 
mIn 

MlN 

MIN 

MlN 
MIN 
MlN 
MlN 
MIN 
MlN 
MIN 
MIN 

m'n 

MIN 
MIN 
MlN 
MIN 

MlN 
MIN 

m*n 

MIN 
MIN 
MlN 
MlN 
MIN 
MlN 
MIN 
MIN 
MlN 
PATH  STATIST-MIN 
MIN 
MIN 
MlN 
MIN 
MIN 
MIN 
MIN 
MIN 
MIN 

mIn 

MlN 
MIN 
MIN 


158 

160 
162 
164 
166 
168 
170 
172 
174 
176 
17B 
180 
182 
1B4 
1B6 
188 
190 
192 
194 
196 
1  ?8 
^00 

202 

204 

?06 
208 
210 

212 

214 
?16 
218 
220 
222 
224 
226 
228 
230 
?32 
234 
236 
238 
240 
2*2 
?44 
?46 
?4R 
?50 
252 
254 
256 
?58 
260 
?6? 
264 
266 
268 
270 
?72 
274 
276 
278 
280 
282 
284 
286 
288 
290 
292 
294 
296 
298 

300 
302 
30* 
306 
308 
310 
312 


- 
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7n. 
71  . 

7?. 

71. 
74. 

7«S. 

7ft. 

77. 


78. 
79. 

80, 


HI  . 
8?. 
BT, 
8*. 
85, 
8*. 


87. 


C 

10S    CONTINUE 

inilT  =  IOUT-l 

TOTAL (KOUNT) =TSUM(NDEST) 

RATIO (KOUNT )=TOTAL< KOUNT) /TOTAL (1 ) 

RVAL*HAT'IO  (KOUNT) 

NUMLK <KOUNT)=IOUT 

IF  (IFTX.EQ.ICOOEO))  CALL  PRTOUT  (7) 

IF  (RATlO(KOUMT) .CjE.?.0. OR. KOUNT. bt.NCHCE)  GO  TO  lib 

c --  — 

C    CALCULATE  ThE  OVERLAPS 

r   __.______------__.---..-----------. 

INUM=IOUl»n.035*n.5 
IF  (INIIM.LT. 1)  INUM  =  1 
KNUMbIOUT»2«INUM 
r ---- 

C    RESET  CENTRAL  Tt  INK  VALUES  TO  VERY  LARGE  NUMBERS 

C — -  —  - 

00    1  in    LL  =  1 »KMUM 
LVAl  =1NIJM  +  LL 
L=ITAH(LVAL) 
1 1  n         TLINK(L)s5000,0 

IF     (IFIX.EQ.IC.ODEO)  )    CALL    PRTOUT     (8) 
GO    TO    «< 
C 

n -_- 

c 
c 
c 

115 


ALL  CHOICES  KFQUESTEU  HAVE  BEEN  FOUND. 

CALL  PRTnUT  TO  DISPLAY  PATH  STATISTICS  AND  PATHS  TO  PREPARE  A  MAP. 


CALL  PPIOUT  (4) 
CALL  PAlHS  (WIDTH) 
RETURN 
END 


MIN 

31* 

MIN 

316 

MIN 

31B 

MIN 

320 

MIN 

322 

MIN 

32* 

MIN 

326 

MIN 

328 

MIN 

330 

MIN 

332 

MIN 

33* 

MIN 

336 

MIN 

338 

MIN 

340 

MlN 

3*2 

MIN 

3** 

MIN 

3*6 

MIN 

348 

MIN 

350 

MIN 

352 

MlN 

354 

MIN 

356 

MIN 

358 

MIN 

360 

MIN 

362 

MIN 

36* 

MIN 

366 

MIN 

368 

MIN 

370 

MIN 

372 

MIN 

374 

MIN 

376 

MlN 

378 

e)        FORSTAT    CROSS-REFERENCE   LISTINGS    FOR  SUBROUTINE    MPATH 

SUBROUTINE  mpath 


ENT  NUMHER 

DEFINITION 

5 

13. 

10 

14. 

15 

?o. 

20 

21  . 

?5 

23. 

30 

26. 

15 

?o. 

40 

3?. 

45 

33. 

SO 

36. 

55 

39. 

60 

43. 

*5 

45. 

7  0 

49. 

75 

53. 

HO 

54. 

85 

56. 

90 

59. 

95 

6?. 

1  0  0 

66. 

105 

70. 

110 

84. 

115 

87. 

REFERENCES 


B6. 

U. 

22. 

58. 

20. 

20. 

2b. 

28. 

23. 

25. 

2b. 

55. 

35. 

3*. 

35. 

A2. 

42. 

**. 

48. 

4V. 

52. 

52. 

5b. 

48. 

69. 

6b. 

6b. 

Bl. 

77. 

28. 


23 

STnTF^ENT  NUMBERS 

DEFINtD 

T 

MAXIMUM  NUMBER  OF 

REFEHtNCES 

3° 

TOTAL    NUMHER  OF 

REFERENCES 

1.3 

AVFRAGE  NL1MRER  OF 

I'FFFKENCES 
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-u 

ul 

JJ 
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tr 

2 
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2 

2 
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f) 


LISTING    OF   SUBROUTINE    PRTOUT 


SUBROUTINE  PRTOUT  (IKOUNT) 

SUBROUTINE  PRTOUT 


PURPOSE  -  TO  PRODUCE  PRINTER  OUTPUT  F OR  PROGRAM  NETSUM  AND  SUB- 
ROUTINE  MPATH. 


THE  GtflRS  SYSTEM 
PURDUE  UNIVERSITY 


A.  KEITH  TURNER, 
CIVIL  ENGINEERING, 
JUNE  1969, 


USAGE—   CALL  PRTOUT  (K) 

WHERE-  K  «  ARGUMENT  USED  IN  COMPUTED  GO  TO  STATEMENT  IN  ORDER 
TO  SELECT  CORRECT  OUTPUT  FOR  EACH  CALL. 


DESCRIPTION  OF  CONTROL  CARDS — 
NONE  REQUIRED 


REQUIRED  SUHROUTINES-- 

NONF 
OPTIONAL  SU9HOUTINF.S-- 

NOME 


INTEGER  SYSOT 

DIMENSION  ITAH(200),  NOOE(20n,7) 

DIMENSION    AREA(3),    FACT0RI4),    NXJ(9800)»    TLINK(9800)»    2(50,50) 

DIMENSION    ZZ(50,50).    TSUM(2500)»    NN(2500) 

DIMENSION  TOTALIVIt  RATIOI7),  NUMLK<7>»  ICODEO) 

EQUIVALENCE  (Z.NIJ),  (ZZ«TSUM) 

COMMON  NXJ,TLINK,Z,ZZ«AHEA,FACTOR»NCOL»NROW,KOUNT»NUMLK»K,N»J«RVA 
I  ,IOUT,VALUE.TnTAL»RATIO,ITAB,NOO£,NHOMElNUEST,NCHCE»IM,IP,IFIX.NF 
?C 

DATA  (ICOOE (I) ,I=1.3)/5HPUNCH.5HmAH   ,5HDEBUG/ 

DATA  IPACK/40000H/»SySOT/6/ 

C — 

C    OUTPUT  ROUTINE  -  OUTPUT  REQUEST  DETtHMINEU  BY  IKOUNT 
C 

GO  TO  <5,ln,lS,25,35,40,SO,60,70,BO,90,V5) ,  IKOUNT 
C    

5      WRITE  (SYSOT, 155) 

WRITE  (SYSOT, 160) 

WRITE  (SYS0T»165) 

WHITE  (SYSOT, 160) 

WHITE  (SYSOT, 170)  AREA 

WRITE  (SYSOT, 175) 

WRITE  (SYSOT. 180) 

WRITE  (SYSOT, 175) 

GO  TO  105 


WRITE  (SYSOT, IR5)  VALUE, FACTOR 

WRITE  (SYSOT, i 75) 

IF  (IP.tU.ICODE (j) )  PUNCH  230,  N,NF AC , VALUE , FACTOR 

GO  TO  105 


w«ite    (Sysot, noi    NCOL.NROW 

00  20  I=1,NR0W 

WRITE  (SYS0T»195) 
WRITE  (SYSOT, 200)  (ZZ(IiJ) tJaliNCUL) 
GO  TO  105 


OUT 
OUT 

-OUT 
OUT 
OUT 
OUT 
OUT 

-OUT 
OUT 
OUT 
OUT 
OUT 
OUT 
OUT 
OUT 
OUT 

-OUT 
OUT 
OUT 
OUT 
OUT 
OUT 

-OUT 
OUT 
OUT 
OUT 
OUT 

-OUT 
OUT 
OUT 
OUT 
OUT 
OUT 
OUT 

-OUT 
OUT 
OUT 
OUT 
OUT 
OUT 
OUT 

LOUT 

AOUT 
OUT 
OUT 
OUT 
OUT 
OUT 
OUT 
OUT 
OUT 
OUT 
OUT 
OUT 
OUT 
OUT 
OUT 
OUT 
OUT 
OUT 
OUT 
OUT 
OUT 
OUT 
OUT 
OUT 
OUT 
OUT 
OUT 
OUT 
OUT 
OUT 
OUT 
OUT 
OUT 


2 

* 
6 

e 

10 
12 

n 

16 
18 
20 
22 
2* 
26 
28 
30 
32 
34 
36 
38 
40 
42 
44 
46 
48 
50 
52 
54 
56 
58 
60 
62 
64 
66 
68 
70 
72 
74 
76 
78 
80 
82 
84 
86 
88 
90 
92 
94 
96 
98 
100 
102 
104 
106 
10B 

no 

112 
114 
116 
US 
120 
122 
124 
126 
128 
130 
132 
134 
136 
138 
140 
142 
144 
146 
148 
150 
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JO. 

k'b 

31. 

3?. 

33. 

34. 

35. 

36. 

37. 

30 

3fl. 

C 

c 

39. 

35 

40. 

41. 

c 

c 

42. 

4fl 

43. 

44. 

45. 

46. 

47. 

4B. 

45 

49. 

C 

c 

5P. 

50 

51. 

5?. 

53. 

54. 

55. 

56. 

57. 

58. 

59. 

55 

6n, 

C 

c 

61. 

60 

6?. 

63. 

64. 

65. 

66. 

67. 

68. 

65 

69. 

C 

C 

70. 

70 

f\> 

l\. 

73. 

74. 

75 

75. 

C 

C 

76. 

80 

77. 

78. 

79. 

eo. 

85 

81. 

C 

c 

8?. 

90 

83. 

C 

c 

84. 

95 

8S. 

Hh. 

ion 

87. 

10S 

WRITE  (SYS0T.P05) 
WRITE  (SYS0T.175) 
WRITE  <SYSOT,210) 
•RITE  (SYS0Ttl75) 
DO  30  I=ltKOUNT 

WRITE  (SYS0T.215)  I . TOTAL ( I ) tRAT 10 ( I ) .NUMLK ( I ) 

WRITE  (SYSOTtl75) 
CONTINUE 
GO  TO  105 


WRITE  (SYSOT.220)  KOUNT 
WRITE  (SYS0Ti225) 
GO  TO  10b 


XI  =  ( FLOAT  <K) /FLOAT (NCOL) >♦ 1.0 

T  X i =X1 

IY1  =  ( (Xl-IXI)OFLOAT  (NCOL) >*0.3 

IF     ( I Yl .NE.n)     GO    TO    45 

1X1=1X1-1 

IYl=NCr>L 

WRITE     (5YSOT.110)     1X1 il Yl iTSUM(K) tK 

GO    TO    105 


WRITE  (5Y5OTil40)  RVAL 
WRITE  (SYS0T.115)  10UT 
WRITE  (SYSOT.120) 
WRITE  (SYSOT.130) 
00  55  IV=l.IOUT 

L=ITAB(lv> 

N=N«J(L) /IPACK 

J=NXJ(L)-N»IHACK 

WRITE  (SYS0Ttl25>  L iNt J, TL INK (L ) 
CONTINUE 
GO  TO  105 


WRITE  (SYSOTil35) 
WRITE  (SYSOTti30) 
DO  65  IV=ltIOUT 

L=ITAB(IV) 

N=NXJ(L)/IPACK 

J=NXJ(L)-N»IPACK 

WRITE  (SYS0T.J25)  L iN. J. TLINK (L ) 
CONTINUE 
GO  TO  10b 


WRITE  <SYSOTtl50>  FACTOR 
00  75  Isj.NROW 

WRITE  (SYS0Til9b) 

WRITE  (SYSnTi200)  (Z ( I t J) . J=l inCoL) 
CONTINUt 
GO  TO  105 


WRITE  (SYS0T.145)  FACTOR 
00  8<S  I«l»NROW 

WRITE  <SYS0Ttl95> 

WRITE  <SYSOTt200)  (Z ( I » J) » J=l iNCOL) 
CONTINUE 
GO  TO  105 


WRITE  (SYSOT.245) 
GO  TO  io5 


PUNCH  ?3b.  NCOLiNROW 

00  jOO  I=j.NROW 

PUNCH  240,  (ZZ(ItJ) ,J«l.NCOL ) 

RETURN 


OUT  152 
OUT  154 
OUT  156 
OUT  158 
OUT  160 
OUT  162 
OUT  164 
OUT  166 
OUT  168 
OUT  170 
OUT  172 
OUT  174 
OUT  j76 
OUT  178 
OUT  180 
OUT  182 
OUT  184 
OUT  186 
OUT  188 
OUT  190 
OUT  192 
OUT  194 
OUT  196 
OUT  198 
OUT  200 
OUT  202 
OUT  204 
OUT  206 
OUT  208 
OUT  210 
OUT  212 
OUT  214 
OUT  216 
OUT  218 
OUT  220 
OUT  222 
OUT  224 
OUT  226 
OUT  228 
OUT  230 
OUT  232 
OUT  234 
OUT  236 
OUT  238 
OUT  24o 
OUT  242 
OUT  244 
OUT  246 
OUT  2*8 
OUT  250 
OUT  252 
OUT  254 
OUT  256 
OUT  258 
OUT  260 
OUT  262 
OUT  264 
OUT  266 
OUT  268 
OUT  270 
OUT  272 
OUT  274 
OUT  276 
OUT  278 
OUT  280 
OUT  282 
OUT  284 
OUT  286 
OUT  288 
OUT  290 
OUT  292 
OUT  294 
OUT  296 
OUT  298 
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C 

C 

88. 

llo 

89. 

115 

90. 

120 

91. 

125 

92. 

130 

93. 

135 

94. 

140 

95. 

145 

96. 

150 

97. 

155 

98, 

160 

99. 

165 

100. 

170 

101. 

175 

102. 

180 

103. 

185 

104. 

190 

105. 

195 

106. 

200 

107. 

205 

108. 

210 

109. 


110. 
111. 

112. 
113. 
114. 
115. 


116. 


FORMATS 


,15) 


215 


220 
225 

230 
235 
240 
245 


(1H  ,5X,I2»1H,.I2.3X,F15.3»20X»I5) 
(1H0.37HTOTAL  NUMBER  OF  LINKS  IN  THIS  PATH  ■ 
(1H0.35HORIGINAL  LINK  VALUES  FOR  THIS  PATH.) 
(1H  ,bX,IS.bX,Ib.bX,I5,5X,U5.3) 
<lHO,8X,lHL.10X,lHN,10X,lHJ.10X,bHTLlNK,//> 
(1H0.34HREVISEO  LINK  VALUES  FOR  THIS  PATH.) 
(1H0.12HTHIS  PATH  IS, F10.3i23H  TIMES  THE  FIRST  CHOICE) 
(1H0,26HSCALED  UTILITY  MATRIX  FOR  ,4A10) 


•  OUT 
OUT 
OUT 
OUT 
OUT 
OUT 
OUT 
OUT 
OUT 
OUT 


FORMAT 
FOHMAT 
FORMAT 
FOHMAT 
FORMAT 
FORMAT 
FORMAT 
FOHMAT 
FORMAT  (1H0,53HDE8UG  PRINTOUTS  REQUESTED,  LISTING  OF    INPUT  OATA  FOOUT 

1R  .4A10»/)  OUT 

FORMAT  (1H1)  OUT 

FOHMAT  (1H  .55X.29H -.-.—.  —  .—_-...-.-..)  OUT 

FORMAT  (1H  ,55X,29HI       THE  GCARS  SYSTEM      I»/1H  »55X.29HI  PUROUT 

10UE  UNIVERSITY  VERSION  I)  OUT 

FORMAT  <///lH  .53HSEARCH  ACTIVITY  -  PREPARATION  OF  THE  UTILITY  SUROUT 

lFACE.,//ilH  ,20X,20HTHIS  ANALYSIS  IS  FOR, 3A10i ///)  OUT 

FOHMAT  (1H  ,51H-— -— -— — — — — —————.--—OUT 

1)  OUT 

FORMAT  (2H  I,8X,lHI»40X.lHI»/illH  1  WEIGHT  I • 12X, 16HFACT0RS  UTILIZOUT 

1ED,12X,1HI,/,2H  It8X,lHl,40X,lHI)  OUT 

FORMAT  (2H  1  ,BX»1H1.40X.1HI,/,2H  I,2X.F5,2,2H  I,4Al0,lHl)  OUT 

FORMAT  (1H1.23HTHIS  UTILITY  MATRIX  OF  »I3»12H  COLUMNS  BY  tI3,46H  ROUT 

10*S  REPRESENTS  THE  CoMBInEO  UTILITY  SURFACE.,///)  OUT 

FOHMAT  (1H0)  OUT 

FOHMAT  (12F10.3)  OUT 

FOHMAT  (1H0,26HTHE  GENERATEO  ALTERNATIVES.//)  OUT 

FOHMAT  (2H  I,8XilHI,12X»lHI,l5XilHI,HX»lHI,/52H  I  CHOICE  I  PATH  TOUT 

10TAL  I  RELATIVE  PATH  I  NUMBER  OF  Ii/illH  I  NUMBER  I.12X.29HI      VOUT 

2ALUE      I    LINKS    I./.2H  I t BX, Ihl , 12X, lHl , 1SX, 1HI » 1 IX, lHl )       OUT 
FORMAT  (2H  I,8X,lHI,12XilHItl5XilHl,llXilHI,/,5H  I    .12, 5H    I  ,FOUT 

110. 3, 6H  I     ,F7.3,4X,1HI,4X,I3,4X,1HI,/,2H  I ,8X, lHl , 12X, lHI , 15X , 10UT 

2HI,11X,1HI)  OUT 

FORMAT  (1H1,27HANALYSIS  FOR  CHOICE  NUMBER  ,13)  OUT 

FORMAT  (1H0,5X,22HMINIMUM  PATH  REOUESTED,/// ,5X,4HN00E, 10X.BHPATH  OUT 

1SUM)  OUT 

FORMAT  (7HNUM8ER  iI2,4H  OF  ,I2,16H  FACTORS, WEIGHT", F5. 2, 4X.4A10 )   OUT 
FORMAT  (27H   COMBINED  UTILITY  SURF  ACE . »43X, 213 ,2H  2)  OUT 

FORMAT  (8F10.3)  OUT 

FORMAT  (1H1.45HSEARCH  ACTIVITY  -  GENERATION  OF  ALTERNATIVES, ,////, OUT 

150H  REPEATED  USE  OF  A  MODIFIED  VERSION  OF  THE  BRITISH, /,*8H  ROAO  ROUT 

2ESEARCH  LABORATORY  MINIMUM  PATH  ALGORITHM. /,41H  GENERATES  AS  MANY  OUT 

3AS  SEVFN  ALTERNATIVES.!//)  OUT 

END  OUT 


300 
302 
304 
306 
308 
310 
312 
314 
316 
318 
320 
322 
324 
326 
328 
330 
332 
334 
336 
338 
340 
342 
344 
346 
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350 
352 
354 
356 
358 
360 
362 
364 
366 
368 
370 
372 
374 
376 
378 
380 
382 
384 
386 
388- 
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h)  LISTING    OF   SUBROUTINE    PATHS 


7. 
3. 

4. 

%. 
6. 

7. 

8. 

<J. 

in. 


li. 
12. 

13. 
1*. 


15. 
16. 
17. 


SUBROUTINE  PATHS  (WIDTH) 

SUBROUTINE  PATHS 


PURPOSE  -- 

TO  OEVELOP  AND  PRINT  A  MAP  SHOWING  THE  LOCATIONS  OF  THE  ALTER- 
NATIVES GENERATEO  BY  SUBROUTINE  MPATh, 

EACH  ALTERNATIVE  IS  OEFINEO  BY  A  STRING  OF  NUMBERS  WHICH  COR- 
RESPOND  TO  THE  CHOICE  NUMBER  (I.E.  FIRST  CHOICE  OEFINEO  BY  STRING 
OF  ONES.  ETC.).   WHERE  TWO  OR  MORE  CHOICES  OVERLAP  OR  CROSSi  THE 
CHARACTER  BELONGING  TO  THE  LOWEST  CHOICE  NUMBER  IS  THE  ONLY  ONE 
PRINTED. 


THE  GCARS  SYSTEM 
PURDUE  UNIVERSITY 


A.  KEITH  TURNER, 
CIVIL  ENGINEERING. 
JUNE  1969. 


MAP 
MAP 
.MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 

•  MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 

•  MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 

-MAP 

MAP 

MAP 

MAP 

NONE  REQUIRED-  ALL  NECESSARY  VALUES  SUPPLIEO  BY  PROGRAM  NETSUM  MAP 

OR  SUBROUTINE  MPATH  VIA  COMMON  STATEMENT.  MAP 

MAP 

._...........--........................................ ............. MAP 

MAP 
MAP 
MAP 
MAP 
MAP 
MAP 

•  MAP 


USAGE  — 


CALL  PATHS  (WIDTH) 

WHERE-  WIDTH  «  WIDTH  OF  MAP  IN  INCHES, 


DESCRIPTION  OF  CONTROL  CARDS-- 


REOUIRED  SUBROUTINES-- 

-NONE 
OPTIONAL  SUBROUTINES— 

-NONE 


DIMENSION  AREA (3) •  FACTOR(A),  SYM(129>.  PSYM<1?9>,  ISYMB(T)        MAP 

DIMENSION  NUMLKI7)  MAP 

DIMENSION  DUM1  (24342),  OUM2I220)  MAP 

DIMENSION  NODE(200,7),  ICODEO)  MAP 

COMMON  SYM, PS YM.DUM1, ARE A, FACTOR, COLS, ROWS, KOUNT,NU«LK,DUM2, NODE, NMAP 

1H0ME,NDEST,NCHCE  MAP 

INTEGER  HOWS, COL S.COl.MAX, COL, PSYM, PS, GRIDtSYM, BLANK, PLUS, MI NUS.XSTMAP 

10P,YSTOH,EYE,SYSOT  MAP 

DATA  BlANK/lH  / ,PLUS/ 1H*/ ,MINUS/1H-/, (ISYMB(I) »!■! |7) /1H1 i 1H2i 1H3.MAP 


11H4,1H5i1H6,1H7/,EYE/1H1/,SYS0T/6/ 
DATA  (ICODE(I) ,I=1,3)/5HPUNCH,5HMAP   .5H0EBUG/ 
WRITE  (SYSOT.90)  NCHCE , COLS, ROWS 

CALCULATE  MAP  WIDTH  IN  CHARACTERS  (COL) 

IF  (WIDTH. LE. 0.0)  WIDTH»10.0 
IF  (WIDTH. GT. 12,7)  WIDTH=12,7 
COL=WIOTH»10, 0*2,0 
PSb (COL-COLS) / (COLS-1 ) 

TEST  GRID  SIZE-- IF  GRID  SMALLER  THAN  2  CHARACTERS  BETWEEN  COLUMNS, 
AN  ERROR  STATFMENT  IS  PRINTEO  AND  PLOT  ABORTED. 

IF  (PS.GT.O)  GO  TO  5 
WRITE  (SYS0T,9S) 
GO  TO  Rb 


MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 


2 

4 

6 

B 

10 

12 

14 

16 

IB 

20 

22 

24 

26 

29 

30 

32 

34 

36 

38 

40 

42 

44 

46 

48 

50 

52 

54 

56 

58 

60 

62 

64 

66 

68 

70 

72 

74 

76 

78 

80 

B2 

64 

B6 

86 

9  0 

92 

94 

96 

98 

100 

102 

104 

106 

108 

110 

112 

114 

116 

118 

120 

122 

124 

126 

128 

130 

132 

134 

136 

138 

1*0 

142 
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18. 
19. 
20. 
21. 
22. 
23. 
24. 


25. 
26. 
27. 
28. 
29. 
30. 
31. 
32. 
33. 
34. 


35. 


36. 
37. 
38. 
39. 


40. 
41. 
42. 
43. 
44. 
45. 
46. 


47. 
48. 
49. 
50. 


51. 
52. 
53. 

54. 
55. 
56. 
57. 
58. 
59. 
60. 
61  . 
62. 


CALCUl.ATF  THE  REQUIRED  CONSTANTS 

COI_MAX=HS«  (COLS-1  >  *COLS 

XST0P=PS*1 

YSTOPaf (PS*1)«6*5)/10 

NLAST»C0LMAX«2 

NLAS=NLAST-1 

GRID»C0LS-1 

TR0W=R0WS*1 

SET  UP  COLUMN  GRID  MARKS 


MAP  144 

MAP  146 


MAP 
MAP 
MAP 
MAP 
MAP 
MAP 
MAP 


C0L=3 

DO  15  J=1.GRID 

DO  10  I=?.XST0P 
SYM(C0L)=8LANK 
10       COL=COL*l 

SYM(COL)oPI.US 
15     C0L=C0L*1 

SYM(1)=HLANK 

SYM(2)=PLUS 

SYM(NLAST)=B|_ANK 

C 

C    PRINT  MARKS  ACROSS  MAP  TOP 

C    ... ............ 

WRITE  (5YSOT.100)  (SYM ( I ) t I»l iNLAST ) 

C    .... .— — .............. 

C    ESTIMATE  VALUES  FOR  ALL  PRINTING  POSITIONS  ON  THE  MAP 

C       THIS  IS  OONE  BY  DEVELOPING  THEN  PRINTING  ONE  LINE  AT  A  TIME 

c    ............................ 

C    LOOP  THROUGH  THE  MAP  LENGTH  .  MATRIX  ROW  BY  MATRIX  ROW 
c    — ................ 

ITOP=l 

DO  80  MM=2iIROW 

IF  (MM.EQ.IROW)  YSTOP«l 

MsMM-1 

c   ___....___.—._...._._. 

c   loop  through  grid,  line  by  line. 

C    IF  KJ  a  1.  LINE  GOF.S  THROUGH  GRID  NODE. 

C    IF  KJ  ■  ?.  SECOND  LINE  IN  GRID  ■  FIRST  ONE  NOT  THROUGH  NODES, 

C    IF  KJ  =  1»4.  ETC.  ADDITIONAL  LINES  SAME  AS  JUST  PREPARED, 

C   .... — ..... 

00  80  KJ=l»YSTOP 

IF  (KJ.GT.2)  GO  TO  80 
DO  20  N=l,129 
20        PSYM(N)=HLANK 
PSYM(2)=EYE 
PSYM(NLAS)=EYE 
IF  (YSTOP.NE.l.AND.ITOP.NE.l)  GO  TO  30 

C ........ 

C   LOOP  ACROSS  MATRIX  COLUMNS 

c .............. 

DO  2S  N=2,NLAS 
25       PSYM(N)=MINUS 
30       C0L*2 

IF  (KJ.E0.2)  GO  TO  55 

C    ... ............ 

C    KJ  »  1,  A  LINE  THROUGH  NODES 

c    ....... ....... 

PSYM(1)»PLUS 
PSYM(NLAST)=PLUS 
DO  50  J»J .GRIO 
IZCOL=COL 
ICOL=COL«XSTOP 
IZ=(M-l)«COLS*J 
!ZP=IZ*1 
00  45  Nol.NCHCE 
KK=NCHCE-N*1 
NLK=NUMLK(KK) *1 
ISW1=0 
ISW2aO 


148 
150 
152 
194 
156 
156 
160 
MAP  162 
MAP  164 
MAP  166 
MAP  168 
MAP  170 
MAP  172 
MAP  174 
MAP  176 
MAP  178 
MAP  180 
MAP  182 
MAP  184 
MAP  186 
MAP  188 
MAP  190 
MAP  192 
MAP  194 
MAP  196 
MAP  198 
MAP  200 
MAP  202 
MAP  204 
MAP  206 
MAP  208 
MAP  210 
MAP  212 
MAP  214 
MAP  216 
MAP  218 
MAP  220 
MAP  222 
MAP  224 
MAP  226 
MAP  228 
MAP  230 
MAP  232 
MAP  234 
MAP  236 
MAP  238 
MAP  240 
MAP  242 
MAP  244 
MAP  246 
MAP  248 
MAP  250 
MAP  252 
MAP  254 
MAP  256 
MAP  258 
MAP  260 
MAP  262 
MAP  264 
MAP  266 
MAP  268 
MAP  270 
MAP  272 
MAP  274 
MAP  276 
MAP  27B 
MAP  280 
MAP  282 
MAP  284 
MAP  286 
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61. 
64, 
65< 
66. 

67. 
68. 
69. 
70. 
71. 
7?. 
73. 
74. 
75. 


76. 
77. 
78. 
79. 
80. 
81  . 
8?. 
8"». 
84. 
8<S. 
86. 
87. 
88. 
89. 
90. 
91. 
9?. 
91. 
9i. 
9^. 
9k. 
97. 
98. 
99. 

ion. 

101  . 
10?. 


103. 


10ft. 

los. 


107. 
108. 
1  ou  , 


00  35  NN=1.NLK 

IF  (IZ.EQ.NODE (NN.KK) >  ISW1«1 
IF  (1ZP.EQ.N0DE(NN»KK)I   ISW2=1 
IF  (ISWl.tO.l)  PSYM(IZCOL)=ISYMB<KK) 
IF  (ISW2.EQ.1)  PSYM(ICOL>=ISYMB(KK) 
IF  (ISWl.EQ.O.OK.ISw2,EQ.O)  GO  TO  *5 
COI.  =  IZC0L 
00  40  I*2,XST0P 

C0L=C0|_»1 
PSYM(COL)=ISYMB(KK) 
CONTINUF 
C0L=IC0L 
GO  TO  75 


KJ  =  ?,  FIRST  INTERMEDIATE  LINE 


55  00    70    J=1,GRID 

?COL=COL*XSTOP 
TZ=(M-1 )*COLS*J 
IZP=IZ*1 
IZOL=IZ*COLS 
IZUR=IZP*C0LS 
00    6b    N=1»NCHCE 
KK=NCHCE-N*1 
NLK  =  NUM(.K(KK)  *1 
ISWl=0 
ISW2=0 
ISw3=o 
ISW4=0 
00    60    NN  =  1 «NLK 

IF     (IZ.EO.NOOE(NN.KK) )     ISW1«1 
IF     (IZP.EQ.NODEINN.KK) )     ISW2=1 
IF     <IZDL.EQ.NODE(NN,KK> )     ISW3M 
60  IF     (IZDR.EQ.NODE (NN.KK) >     ISW*=1 

ITEST=ISW1 
ITFST=ISW2 
ITFST  =  ISn*3 
ITFST=ISW4 

IF     (ISW) .EQ.l. AND.ISW3.EQ.il     PSYM(COL) 
IF     (ISW2.E0.1,AN0.ISW*.E0.1)     PSYMdCOL 
65  CONTINUF 

70  COL=ICOL 

75  ITOP=0 

r --- — - — - — - 

C    PRINT  OUT  LINE  JUST  PREPARED 

r   ------- — --.-_......--_..-.-.._-.___ 

80     WRITE  (SYSOT.100)   (PSYM(I) tlel.NUAST) 

r 

C    «AP  FINISHED,  PRINT  MARKS  ACROSS  MAP  BOTTOM 

f . 

WRITE     (SYSOT.100)     < SYM < I ) , I =1 ,NLAST ) 
85             RETURN 
C         FORMATS       


40 


4S 

loo 


FORMAT  (1H0.//.31H  MAP  OF  GENERATtO  ALTERNATIVE 
10ICES  ARE  SHOWN. ./.45M  EACH  IS  LOLATFD  BY  A  STR 
?H,/,33H  CORRESPOND  TO  THE  CHOICE  NUMBER. ,//,42H 
^"E  DEFINED  ON  A  NETWORK, /,8H  HAVING  .I2,12H  COL 
4..///) 

FORMAT  (25H  TOO  MANY  COLUMNS  IN  GHIOl 

FORMAT  (1H  ,1?9A1) 

END 


MAP  ?BB 

MAP  290 

MAP  292 

MAP  29* 

MAP  296 

MAP  298 

MAP  300 

MAP  302 

MAP  30* 

MAP  306 

MAP  308 

MAP  310 

MAP  312 

MAP  31* 

MAP  316 

MAP  318 

MAP  320 

MAP  322 

MAP  32* 

MAP  326 

MAP  328 

MAP  330 

MAP  332 

MAP  33* 

MAP  336 

MAP  338 

MAP  3*0 

MAP  3*2 

MAP  3** 

MAP  3*6 

MAP  3*8 

MAP  350 

MAP  352 

MAP  35* 

MAP  356 

MAP  358 

MAP  360 

MAP  362 

MAP  36* 

■ISYMB(KK) 

MAP  366 

(•ISYMB(KK) 

MAP  368 

MAP  370 

MAP  372 

MAP  37* 

MAP  376 

MAP  378 

MAP  380 

MAP  382 

MAP  384 

MAP  386 

MAP  3BR 

MAP  390 

MAP  392 

MAP  396 

S.»//,1H  «I?,19H  CHMAP  398 

ING  OF  NUMBERS 

WHICMAP  400 

THE  ALTERNATIVES  AMAP  *02 

UMNS  BY  ,I2,6H 

ROWSMAP  404 

MAP  406 

MAP  *08 

MAP  *10 

MAP  *12 
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i)  FORSTAT    CROSS-REFERENCE    LISTINGS    FOR    SUBROUTINE    PATHS 

SUBROUTI"t  PATHS 


STATFmEM  NUMBER    DEFINITION 


REFERENCES 


5 

IK. 

lb. 

]  P 

29. 

27. 

lb 

31. 

26. 

20 

43. 

42. 

?S 

4R. 

47. 

30 

49. 

4b. 

15 

6?. 

63. 

4  0 

7?. 

70. 

45 

7?. 

SB. 

50 

74. 

53. 

SS 

76. 

50. 

60 

93. 

89. 

f.b 

100, 

82. 

70 

101. 

76. 

7b 

102. 

7b. 

HO 

103. 

37. 

MS 

10*. 

17. 

90 

106, 

FORMAT 

10. 

>35 

107. 

FORMAT 

lb. 

100 

108, 

FORMAT 

3b. 

68. 


40. 


41. 


103.   104. 


20 

STATEMENT  NUMBERS 

DEFINED 

3 

MAMMnM  NUMHER  OF 

REFERENCES 

?s 

TOTAL    NUMBER  OF 

REFERENCES 

1  .3 

AVI-  ■  A'iF-  NUMRER  OF 

REFERENCES 
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1/1 
aJ 

s 

o 

U 

o 

2 

2 

-z 

UJ 

LkJ 

X 

CE 

j: 

u. 

a. 

UJ 

UJ 

UJ 

IX 

a 

a 

u. 

u. 

u 

o 

c 

III 

X 

^ 

3. 

aJ 

a 

'1 

i 

2 

* 

3 

3 

_ 

•^ 

2 

* 

J 

T 

u 

— 

C 

1 

_p 

< 

1 

a 

t 

1/1 

£ 

>- 

<J 

t\ 

p. 

(*> 

r. 

if 

f\. 

f^ 

-a 

»-,         (V    s 


«  z  Ttn^-(v      »— c 

UJ  <  _J  jjxjujuh 

a  jo  oodu><ii 

<  X  O  U  U  Q  Q  It'  II  15 


uj  rr  »-  _j  h-  w  uikui      h  *  t-iaa. 
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j)  EXAMPLE     OF    CONTROL    CARDS    DATA    FOR    PROGRAM    NETSUM 


2    12    24    32 

2    512.7 

Cr 

IMOIIHS       PUNCH 

1.0       1.0 
1001 

EARTHWORK    COST    FACTOR. 

NORTH 

25    33 

1 

13.368 

13.130 

3. 

U71 

-4.207 

-10.072 

-9.4M3 

-11. 

249 

-1  ". 

725 

-25.202 

-20.409 

-21. 

13') 

-4.460 

-11.345 

-8.477 

283 

10, 

,574 

9. 755 

22.193 

32, 

168 

27.301 

-50.099 

-  35.816 

-40. 

,023 

-24, 

,723 

33.917 

24.74ft 

18.974 

12, 

.  164 

-4.7H3 

-9.860 

-5.301 

-9, 

,436 

-19, 

,185 

-18. 30ft 

-11 .214 

-11. 

S09 

-6.624 

-2.695 

-2.047 

7, 

,7]1 

1  <♦  , 

,338 

33.669 

32.947 

38, 

,  J30 

33.321 

24.661 

23.658 

-26, 

,647 

19, 

,3111 

-1.712 

26.445 

31. ^67 

30, 

121 

-l.lft3 

-7.9hl 

-.223 

-6, 

,  799 

-9, 

,033 

-1.111 

2.0<tb 

5, 

,626 

1 0.067 

10.322 

.120 

29, 

,463 

36, 

,489 

42.922 

30.004 

34, 

,715 

3  0.2)1 

-23.000 

-36.05) 

-29, 

,639 

15, 

,093 

-55.166 

15.450 

18.932 

s, 

,05ft 

-5.100 

-7.570 

-1.612 

5, 

,563 

-1  , 

,035 

-.256 

10.751 

V, 

,915 

19.624 

19.459 

17.355 

39, 

,304 

33, 

,689 

33.591 

19.361 

30, 

,891 

26.754 

.051 

5.128 

-37, 

,740 

-58, 

,459 

-64.740 

3.257 

7.293 

17, 

,4U 

-6.597 

5.165 

-3.256 

,137 

5, 

,170 

3.989 

12.820 

-2, 

,392 

9,004 

14.054 

3  3.24  0 

37, 

,194 

35, 

.095 

-3. 101 

3.734 

^, 

,883 

12.076 

6.078 

-R2.423 

-57, 

,055 

-30, 

,245 

21.597 

-8.275 

-4.875 

-10. 

"91 

.824 

-19.419 

20.608 

7, 

,701 

,501 

7.077 

15.218 

^i, 

,548 

39.537 

28.989 

31.397 

-5, 

,633 

-B, 

,135 

-10.090 

10.669 

-47, 

,546 

-26.037 

-65.027 

?6,061 

56, 

,332 

*5, 

,727 

23.392 

-21.775 

-16.293 

-8, 

,688 

-6.275 

-15.369 

-10.261 

-6, 

,147 

3, 

,427 

8.672 

1.964 

16, 

,626 

29.475 

-11.181 

-8, 144 

21, 

,474 

-27, 

,340 

-38.731 

-44.246 

-2«. 

,  144 

-69,558 

-61.141 

-44.797 

57, 

,025 

38, 

,419 

-18.644 

-33.557 

-29.221 

-21. 

,666 

-25.050 

-27.246 

8.59ft 

-9, 

,098 

4, 

,247 

-5.103 

-15.065 

?2, 

,773 

-2.190 

-1.234 

■36.670 

-22, 

,639 

-49, 

,48ft 

-27.142 

-27.987 

-74, 

,262 

-48.969 

34.499 

49.9J8 

53, 

,873 

»n 

,954 

29.560 

-31.500 

-28.825 

-22, 

,631 

-16.594 

3.555 

1.553 

-9, 

,176 

-15, 

,752 

-11.154 

12.331 

20, 

,842 

33.719 

27,998 

-14.537 

-48, 

,  740 

-29, 

,  144 

-39.529 

-65.788 

,804 

-27.176 

29.709 

40.404 

15, 

,766 

29, 

,246 

16.206 

-21.911 

-17.129 

-18, 

,050 

-11.453 

-7.488 

-4.406 

-22, 

,138 

1, 

,825 

26.143 

43.319 

40, 

,795 

48.3?4 

-26.833 

-76.765 

-66, 

,567 

-74, 

,803 

-69.303 

-49.883 

21  , 

,271 

35.903 

39.411 

14.814 

39, 

,785 

29, 

,393 

21.203 

-14.791 

-6.557 

-11, 

,900 

-25.912 

-27.877 

-30.396 

-, 

,049 

29, 

,057 

36.630 

48.915 

S3, 

,870 

34.531 

-81,ft44 

-78.096 

-80! 

,261 

-23, 

,078 

-50.929 

24. 740 

26, 

,065 

41.81? 

48.005 

43.640 

39, 

,989 

35, 

,131 

19.068 

-7.603 

-ft. 231 

-37, 

,197 

-33.856 

-22.424 

3.513 

15, 

,964 

43, 

,777 

52.459 

61.180 

108, 

,396 

-107.051 

-92.S99 

-88.317 

-86, 

,421 

-21, 

,388 

8.555 

1ft. 281 

25, 

,127 

48,155 

45,561 

41,342 

34, 

,612 

26, 

,919 

16.508 

-13.532 

-7.573 

-28, 

,391 

9.116 

17.679 

28.818 

18, 

,616 

59, 

,023 

74.861 

85.030 

30, 

,49ft 

-fl6,9}4 

-89,743 

-87.156 

-75, 

,454 

-19, 

,024 

2.  789 

1ft. 263 

33, 

,  iHb 

47.208 

42,163 

33.8ft9 

17, 

,004 

18, 

,070 

15.720 

-41.168 

-.98* 

15, 

,0  28 

31.430 

41.259 

39.595 

52, 

,388 

67, 

,487 

89.543 

92.729 

-P6, 

,751 

-"0.728 

-80.941 

-M.  163 

-51, 

,995 

-47, 

,172 

10.960 

12.872 

46, 

,624 

43.98f| 

39.540 

31 .438 

11, 

,668 

9, 

,018 

-10.373 

29.094 

34.071 

38, 

,6*9 

36.165 

33,?41 

■57,  7?8 

60, 

,397 

63, 

,024 

59.n5l 

92.977 

-78, 

,9n9 

-84.6^8 

-73.262 

-49.602 

-25, 

,225 

-43, 

,697 

38.672 

41 .0«0 

31, 

,^07 

19.246 

7.?P7 

-'-7. 9^9 

-52, 

,139 

-36, 

,640 

-.502 
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isi.  n3i 

44.280 

43.301 

4  4.140 

29.191 

4  6, 

,  706 

62, 

6n« 

47, 

138 

51.  135 

86.653 

-/3.095 

-89.798 

-24.427 

2. 

,264 

33. 

642 

2. 

666 

-49.104 

-10.012 

-49.206 

-70.941 

-66.630 

4, 

,574 

7. 

8?2 

-41  , 

4  76 

8.327 

68.022 

-8.240 

5H.094 

43.686 

3  3.470 

46. 

,705 

18. 

,091 

7, 

352 

2.874 

-10. MO 

-87.266 

-39,28? 

-3.526 

■30, 

,741 

56, 

,041 

66, 

,231 

62.581 

-51.801 

-71,b33 

13.196 

9.787 

-47, 

,906 

-47, 

,862 

8, 

,317 

10.732 

■51.681 

67,384 

30.071 

-9,180 

28,835 

-7, 

,537 

-, 

,403 

-3, 

,125 

-12.929 

-AS. 155 

-51.204 

-20,869 

46.061 

79, 

,696 

66, 

,853 

65, 

,506 

57.835 

36.075 

-27.555 

11.628 

8.433 

-, 

,216 

-41, 

,397 

-39. 

,144 

-45.126 

-34.272 

-59.298 

-33.579 

-47,004 

-24.636 

-15, 

,515 

-19, 

,  790 

-19, 

,963 

-37.521 

-80.202 

-19.112 

45,220 

77.3)8 

6*. 

,335 

55, 

,826 

«4, 

,456 

44.984 

33.271 

31.308 

27,656 

13.861 

11, 

,359 

-42, 

,843 

-38, 

,836 

-10.715 

8.007 

-60.214 

-86.588 

-94.789 

-89.?47 

-98, 

,118 

-46, 

,731 

-58, 

,826 

-78.788 

-"6.934 

38.552 

46.296 

49.976 

45, 

,612 

50, 

,371 

4", 

,852 

34.394 

29,572 

18,098 

13,011 

20.941 

15, 

,035 

-47, 

,111 

-17, 

,832 

-38.39) 

31.107 

30.548 

-7,890 

-19,402 

-77.?92 

-87, 

,768 

-81, 

,840 

-78, 

,302 

-7,254 

41.612 

-30.372 

16.634 

36.517 

40, 

,399 

35, 

,699 

20, 

,998 

19.684 

15.661 

16.555 

21.627 

27.115 

20, 

,569 

-39, 

,984 

1, 

,842 

1 0. 1  IB 

40.924 

30.614 

31,406 

49,00? 

-74.Q46 

-B7, 

.783 

-1, 

,783 

4, 

,728 

18.89? 

30.499 

25.648 

31.938 

30.983 

?fl, 

,016 

24, 

,530 

16, 

,097 

16.519 

18,990 

25.173 

18.234 

12.835 

9, 

,408 

-15, 

,977 

-13, 

,171 

22.953 

-51.703 

31.029 

24.190 

16.193 

-2.539 

-57, 

,629 

-27, 

,753 

15, 

,221 

18.876 

32.169 

31.319 

27.213 

23.725 

11, 

,9()0 

15, 

,611 

9, 

,777 

4.516 

4.273 

1.770 

5.224 

12.425 

13, 

,361 

9, 

,865 

-24, 

,787 

40.200 

■389.918 

42.909 

?2.765 

14.219 

6,918 

6, 

,89? 

-50, 

,856 

-54, 

,997 

-18.599 

8.345 

17.402 

22.379 

8.199 

3 , 

.953 

7, 

,7)5 

1, 

,192 

9.097 

1.918 

3.368 

2.104 

6.030 

8, 

,426 

12, 

,984 

-11, 

,060 

-19.482 

52.748 

33.258 

18.811 

3.767 

-1.319 

,689 

-10, 

,819 

-23, 

,738 

-43.441 

-32.226 

-2.521 

2.4?6 

-7.603 

-2, 

,133 

,316 

,828 

-3.911 

-4.438 

21.438 

-1.119 

7.151 

3, 

,460 

15^ 

,158 

1, 

,465 

-4.352 

43.879 

51.375 

31,128 

12.112 

2.?56 

-1  , 

,346 

-8, 

,547 

-  i, 

,845 

4.129 

-5.42b 

-16.417 

-1.846 

24.667 

|7, 

,849 

2, 

,609 

-  i, 

,067 

-8.386 

-14.665 

-10.645 

-1 0,764 

-11.983 

4  , 

,716 

3, 

,566 

3, 

,446 

-2.413 

58.431 

34.007 

35,945 

?6,683 

39.059 

?4( 

,360 

3, 

,426 

-33, 

,840 

4.996 

6.491 

-22.334 

-?0.9o9 

-6.099 

30. 

.371 

32, 

,751 

24, 

,16? 

14.677 

-16.127 

-19.647 

-1«.231 

-22.?77 

-?0, 

.150 

,947 

-11, 

,411 

14.564 

61.233 

82.668 

32,"»12 

45.215 

36,564 

24, 

.167 

-3, 

.063 

,058 

-12.660 

3.867 

.404 

-3.385 

-15,?45 

-39, 

,840 

-16, 

,693 

-18 

,238 

-22.334 

-?B,e>25 

-28,221 

-27,09] 

-8.281 

-, 

,325 

5, 

,462 

9, 

,56? 

18.113 

39.616 

21.248 

21.280 

29.995 

34.409 

34, 

,296 

)6, 

,914 

-7, 

,243 

-15.180 

-1 7.255 

-18.323 

-12.143 

-22.305 

-27, 

,670 

-38, 

,251 

-31, 

,155 

-27.34R 

-25.522 

-22.898 

-19.813 

-3.522 

1  T, 

,431 

19, 

,629 

39, 

,105 

38.786 

26.750 

-10.137 

.591 

14.1?6 

?6.799 

39, 

,073 

22, 

,600 

23, 

,350 

12.290 

-24.165 

-23.607 

-23.81? 

-19.703 

-17, 

,99b 

-28, 

,969 

-11, 

,484 

-20.000 

-19.497 

-15.501 

-6,38? 

2.926 

1, 

,470 

22, 

,816 

33, 

,065 

33.807 

-14.399 

-3.812 

-16.440 

-6,846 

7.545 

1«, 

,374 

2), 

.679 

24, 

,337 

10.580 

-15.044 

.067 

2.1?4 

14.035 

14, 

,520 

-25, 

,474 

-30, 

,68  6 

-8.b47 

-.639 

-.174 

5,523 

4. 740 

?, 

,578 

5, 

,839 

)2, 

,09? 

-11.603 

-32.714 

-33,552 

-26.327 

-17.395 

-16.759 

9 

.812 

20 

.577 

12, 

.111 

-.306 

-17.5?2 

25./ 16 

34,329 

34.126 

40 

.  127 

23 

.902 

-31 

,410 

4.619 

-1 .089 

-7,  121 

-5.5no 

-1.462 

-  v 

,9  tH 

,9?4 

-5, 

,266 

-1 8.730 
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-59.458 

-55.122 

-43.643 

33.636 

-?3.?47 

-9.275 

1.742 

5.79? 

-11.496 

5.597 

16.662 

43.287 

52.00b 

51.043 

41.253 

11.633 

-?4.?55 

.002 

.3*6 

1.7S0 

1.098 

-,39H 

-4.433 

-9.435 

-3.086 

PAVEMFNT    CONSTRUCTION 

COST    FACTOR 

NORTH 

?5 

33    2 

5.283 

5.281 

b.288 

5.357 

5.3H5 

5,697 

5.358 

5.300 

■5.474 

5.474 

5.467 

5.7?8 

5.S07 

5.707 

6.9S8 

5.770 

5.300 

5.481 

4.892 

2.413 

5.900 

1.965 

1.300 

1.89? 

1.747 

5.219 

5.224 

5.289 

5.38? 

5.749 

5.772 

5.300 

5.26) 

5.545 

5.611 

5.576 

5.736 

5.615 

7.199 

7.411 

5.761 

8.000 

5.288 

4.496 

2.537 

2.?52 

1.919 

1.300 

1.819 

2.003 

5.200 

5.276 

5.200 

5,4*8 

5,900 

5.300 

5.900 

5.273 

5.71B 

5.300 

5.540 

5.880 

5.939 

8.000 

5.732 

5.674 

5.300 

4.861 

l.t>00 

3.288 

7.?36 

2.802 

4.661 

2.282 

2.325 

5.200 

5.300 

7.179 

5.715 

7.188 

7.023 

5.388 

5.704 

8.636 

8.494 

5./15 

5.440 

5,86b 

5.532 

5.4H5 

5.300 

4.532 

4.471 

2.799 

3.172 

7.186 

5.212 

4.880 

5.203 

5.331 

7.280 

7.441 

7.374 

7.283 

5.64  6 

6.949 

5.798 

6.397 

8.673 

6.18b 

5.300 

5.7?8 

5.759 

5.3*5 

5.158 

4.558 

2.221 

2.79b 

3,084 

3.214 

2.641 

1.500 

4.876 

5.235 

5.331 

5.737 

5. 860 

5.571 

5.379 

5.?14 

5.200 

6.9?9 

7.407 

5.804 

5.707 

5./00 

5.300 

5.772 

5.300 

5.434 

5.030 

1.985 

3.036 

2.162 

6.03? 

5.^47 

5.355 

5.015 

5.300 

5.151 

5.606 

5.439 

5.284 

5.277 

5.?13 

5.27? 

5.162 

5.245 

5.328 

5.300 

5.420 

5.426 

5.434 

5.900 

5.466 

5.300 

7.094 

8.000 

2.297 

2.590 

5.?21 

1.500 

5.900 

1.300 

5.073 

5.625 

5.656 

5.691 

5.456 

5.219 

5.218 

4.856 

2.551 

5.  149 

5.900 

5.300 

5.7?8 

5.404 

5.09b 

5.146 

7.553 

3.261 

2.717 

2.701 

4.8?0 

2.504 

4.156 

1.868 

4.640 

4. 825 

5.807 

5.960 

7.486 

5.SHO 

5.481 

5.300 

4.705 

5.200 

5.406 

5.728 

5.689 

5.376 

5.092 

4.726 

2.9?2 

1.300 

2.070 

4.971 

4.945 

1.300 

5.300 

5.092 

5.300 

5,403 

5.109 

5.758 

6.006 

7.672 

7.476 

5.903 

5.406 

2.613 

2.796 

5.300 

5.279 

5.316 

4.7?6 

1.300 

4.613 

2.170 

?.8H5 

5.900 

4.821 

2.390 

5.300 

5.900 

5.3  7? 

5.714 

5.364 

5.332 

5.495 

7.122 

7.437 

6.964 

5.?51 

5. 154 

5.300 

4.99? 

5.294 

5.295 

5.037 

4.7M1 

5.900 

3.174 

5.900 

?.465 

2.467 

2.125 

1.500 

4.399 

5.133 

5.383 

5.348 

5.34? 

5.338 

5.345 

5.355 

5.253 

4.760 

1.300 

2.958 

5.604 

5.338 

5.300 

5.300 

5.055 

4.973 

5.?17 

5.286 

4.991 

2.281 

1.869 

2.108 

2.318 

2.740 

5.125 

5.300 

5.300 

5.300 

5.300 

5.079 

5.200 

5.083 

1.300 

4.684 

5.300 

8.184 

5.835 

5.300 

5.342 

5.021 

5.047 

5.035 

4.9?1 

5.036 

4.494 

2.034 

2.016 

5.300 

4.  788 

5.330 

5.338 

5.300 

5.328 

5.328 

5.161 

1.300 

5.575 

5.300 

5.059 

5.693 

8.495 

5.884 

5.867 

5.343 

5.300 

4.839 

3.004 

3.054 

4,863 

2.381 

1.993 

2.268 

4.628 

4.551 

5.110 

4.662 

4.981 

5.123 

5.157 

2.388 

4.820 

7.371 

5.300 

5.668 

5.200 

5.300 

8.  104 

5.300 

4.969 

5.202 

3.000 

1.300 

4,619 

4.576 

1.500 

5.300 

4.434 

4.065 

2.512 

2. 555 

5.900 

1.300 

5.064 

5.551 

5.060 

5.610 

6.098 

7.248 

5.839 

5.279 

5.603 

5.4?0 

5.081 

5.082 

4.944 

1.300 

2.087 

?.03O 

2.0*3 

4.388 

?.?93 

2.172 

4.28? 

2.041 

5.900 

2.501 

5.900 

5.1?B 

5.  199 

B-52 


5.425 

5.  /25 

5.372 

S.47? 

5,f-83 

5,144 

5.  3«»7 

5.900 

2.615 

2.229 

2.299 

2.247 

1  .845 

2.233 

4.7(14 

5.097 

4.877 

2.114 

5.900 

1.838 

2.317 

2.331 

5.19b 

5.431 

5.418 

5.300 

5.413 

5.133 

5.9no 

4.696 

5.20  0 

5.285 

2.360 

2.1*8 

4.453 

2.339 

2.076 

4.428 

4.938 

5.300 

5.569 

S.501 

4.546 

1.300 

4.0*6 

1.500 

2.375 

5.068 

5.426 

5.386 

5.900 

2.423 

3.J15 

3.000 

2.278 

2.171 

1.972 

1  .300 

1.969 

2.604 

4.955 

4.720 

4.966 

5.413 

5.553 

8.000 

5.771 

5.300 

2.388 

5.148 

5.426 

5.301 

5.900 

5.202 

5.773 

5.900 

4.98b 

3.033 

5.900 

3.368 

2.831 

2.118 

1.300 

5.900 

2.601 

2.886 

5,161 

5.352 

5.300 

5.519 

7.357 

5.692 

8.000 

5.478 

7.200 

7.612 

6.634 

2.359 

5.030 

5.900 

2.352 

2.471 

5.026 

2.918 

2.536 

2.976 

4.738 

5.043 

5.110 

4.864 

5.183 

5.900 

5.577 

5.268 

5.346 

7.098 

8.000 

7.660 

7.b45 

6.362 

5.509 

2.528 

1.300 

4.816 

5.169 

2.174 

2.214 

2.350 

4.763 

2.167 

2.260 

1.500 

2.487 

2.584 

5.034 

7.232 

7.336 

7,300 

5.355 

3.273 

5.705 

8.000 

7.752 

7.458 

6.775 

6.306 

1.500 

1.679 

2.637 

4.982 

1.837 

2.540 

2.162 

1.908 

1.909 

4.512 

4.516 

2.731 

4.930 

4.626 

7.014 

7.046 

6.877 

5.132 

3.035 

7.141 

7.503 

6.204 

5.483 

4.931 

2.614 

3.817 

1.500 

2.268 

4.633 

2.517 

2.561 

1.598 

1  .732 

1.755 

2.014 

2.379 

6.827 

8.000 

2.576 

4.9Q2 

5.166 

5.300 

2.941 

2.650 

5.531 

7.127 

5.634 

3.411 

6.436 

4.599 

1.880 

1.990 

1.300 

4.967 

1.831 

1.710 

1.704 

2.232 

1.500 

2.382 

2.175 

2.275 

5.012 

2.445 

8.050 

6.785 

7.615 

7.795 

2.969 

3.141 

3.384 

5.293 

5.273 

6.757 

3.153 

2.206 

2.124 

2.447 

5.324 

3.053 

2.942 

2.232 

2.400 

4,894 

6.4)8 

2.285 

3.000 

1.500 

1.300 

7.393 

3.365 

5.393 

3.210 

8.000 

3.159 

6.619 

6.756 

3.039 

8.000 

5.263 

4.572 

2.464 

1.500 

6.641 

1.50O 

4.776 

5.0  09 

6.248 

3.937 

8.000 

3.000 

4.906 

4.757 

4.807 

5,500 

4.625 

4.838 

5.219 

5.260 

3.250 

5.20O 

3.098 

3.200 

7.203 

7.343 

1.500 

4.332 

5.371 

5.184 

5.464 

5.515 

5.352 

1.500 

5.087 

5.200 

6.615 

4.59) 

5.900 

4.904 

4,730 

2.179 

l.^OO 

5.300 

5.200 

4.9*9 

1.500 

7.056 

7.212 

7.502 

7,342 

4.894 

2.215 

2.306 

2.275 

5.422 

5.597 

5.715 

6.800 

7.171 

5.797 

7.082 

1  .500 

2.401 

2.579 

6.213 

2.570 

2.708 

6,923 

8.000 

8,000 

2.519 

6.660 

1.500 

5.200 

5.784 

5.0  74 

4.686 

4.878 

4.634 

5.  ?64 

5.433 

8.U00 

5.79? 

5.560 

5.460 

5.338 

4.995 

4.456 

5.137 

6.635 

8.000 

1.500 

6.699 

5.489 

5.300 

4.473 

5.200 

5.101 

5.147 

5.216 

5.200 

5.2P5 

5.127 

5.108 

5.737 

5.847 

5.642 

5.530 

5.440 

5.421 

5.285 

5.286 

5.055 

5.300 

5.945 

6.199 

6.793 

5.096 

5.900 

2.614 

4.665 

5.068 

5.225 

5.277 

5.279 

5.29? 

5.146 

5.131 

5.174 

5.392 

5.507 

5.492 

5.433 

5.?94 

5.2M3 

5.268 

5.351 

5.327 

5.900 

7.396 

5.618 

5.200 

5.323 

5.502 

3.042 

5.300 

5.060 

5.285 

5.288 

5.294 

5.046 

4.973 

2.959 

5.501 

5.473 

5.431 

5.333 

5.419 

5.100 

5.2o3 

5.255 

5,720 

5.«08 

5.200 

5.594 

5,546 

5.402 

5.286 

5.173 

5.4M7 

5.056 

5.200 

5.277 

5.300 

5.322 

5.131 

5.300 

5.186 

5. 084 

k) 
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EXAMPLE    OUTPUT 


I        THE  GCAHS  SYSTEM      I 
I  PURDUE  UNIVERSITY  VERSION  I 


SEARCH  ACTIVITY  -  PRE Para T T UN  OF  THE  UTILITY  SURFACE, 
THIS  ANALYSIS  IS  FORnORTm 


WEIGHT    I 
I 


H  fCTnRS    IJTIL  IZEO 


1.01  I  EARTHWORK  COST  FACTOR 


l.On  I  PAVEMf NT  CONSTRUCTION  COST  FACTOR 


CONTOUR  MAP  OF  THE  constructed  UTILITY  SUHFACE 

MAP  IS  DEVELOPED  FROM   «?5  GRIODED  VALUES. CONS  1ST InG  OF   25  COLUMNS  BY   33  ROWS. 


CONTOUR  LEVELS  CALCULAtED 

GRIOUEO  OATA 

LINEAR  INTERPOLATION  wtTHlN  GRID  SQUARES 

OBSERVED  VALUES 

ROWS  =   33    COLS  =   25 

ZRAR  =       5.370    /STG  = 


VALUES  EMPLOYED 
ZMIN  =        .(125 
GRIO  =      5    XINC 
CON  =     10    PS  = 


2.2BR    WAXZ  = 


11.651    CI 
.200    YINC  = 
COLMAX  =    121 


111651    MINZ 


.025 


1.163 
.333 


CONTOUR  SYMBOL  TABLE  IS  AS  FOLLOWS 

(ONLY  EVEN  SYMHOI.S  PRINTED  IN  BUOY  OF  MAP) 


SYMBOL 


VALliE 


NO  DATA 
BELOW   LOw 


0 

.0252 

TO 

1.1877 

1 

1.1H77 

TO 

2.3503 

2 

?. 35(i3 

TO 

3.5128 

3 

3.512B 

TO 

4.6754 

A 

4.675* 

TO 

5.3379 

5 

5.B37V 

TO 

7,0004 

6 

7  .  O  0  O  4 

TO 

B.1630 

7 

B.1630 

TO 

9, 3255 
10.*BP1 

B 

9.3255 

TO 

V 

10.4BR1 

TO 

11.6506 

♦ 

11  .65"6 

OS       OVr'R 
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CONTOUH  MAP  OF  iME  CONSTRUCTED  UTILIT  SUHFACE 

MAP  IS  DEVELOPED  H-Of   B2J  bHJUDED  VALUES .CONSIST ING  OF   25  COLUMNS  HT   33  R0»S. 

.444444444444444»44444'.S'^SbbSb444bb5bbbSbbSbbbbbbbbbbbb6b5bSbSSbbS5566O66'>5b5'>5444b5555554443321234565432llll000111111111 
4444444444444444444  4444  6666666         6  444   2  2   444  72     00000        1 

5   4444444, 4, ,4, ,4  44444  6666b6  666    666  666     4444   22222   22?     0000OO0     01 

•  5  444444444444a  444444  6666  666  66  M  6  4444  222  00O00000  00000 
5  4444444  4*4  666  666666  666  66  444  222222  22?  0 
5             44444         44«         444         4444                66     66666     66666    44   2222222   44  44  222222  222       0 

.5              4444       44444       4444     44444444             666    666                444  2  0   222  4  6  6  4  2?   44  2?    22? 

5  4  6  44  666666  666  6666  4444  2?  2222  4  6  6  44  4444  ?  2 
44         666        6o*6666              66       666                666666            4444444   22  22222  4  6  6   4  444444444444442 

•4444     66666      66  hMtat    144     6   BB8BB8  b                                                                                     44444444    2222222  4  666    444444  66 

5      6666  6666   666666666*6   '.      66   BBBB    66                               444           22222222   44  444    444  6 

666666666    6666666    666666          66  BOB    666   44                         444   ???2?2222222?2?2222      222   44  6 

«7                 66      C6666         h6   BBH   6     4444                       444   ?          ?2?2222????          ??  »  5> 

66666     666                               4444   ?         2222222            ?22   44  5 

66666b66666666                                44444   ?       2222222   44444444444  44*  5 

666666   666                                   444444   2       ?2?2222  4  6666666  4 

6666666                                      4444444       44   22  2   4      6   44  44  44444444 
444                                          666  4  ?  22   444   6  4   2  44  6  4   22   44 
1444444444444444444444                                66    BBB6   4   2222   44   44  2?464   02   4 

444444     4,444444                 4               666666666666  4   222     44444   ?2?   4  4  2??   44 

44444   72?   4444                   444             66  66   4444444    2    4  4  4  4  444 

4444   22  2?   44                  44444         4  6   B  6   ????222?2    444  4   222    4   2?  ?  4444444 

4444    2?    444                    44444444444444     4   22?                       4444    2  44444444 
444444     44444                     444444444         22222?   444444    2??22   4444444444  4444    444 

4444444444444                        444444    222   00     2  44    444  22  0  ?  4            4  44 

4          4444                   444    44444    2222    ??    44    44    22222  44  4 

22222?   44                  4   222   4  44   22222?   44       44   ??     44  » 

2       24                  42024  44   222    44  6666    4   2  ?   4     666          6  5 

4   ?2??2?   44                   444   2   44444444      2   44   444   22   22  4  5 

44         444                   4444444444444    44   ??         222?     2  44  5 

444444444                     44   66  6   4  4   6  6   2   ?2?2?2        00  2   444  5 

444444                              666666   4     66  4  2                      22   444  5 
44                    66666666666666666    666b  4  2 
44                   6666666666666666666666666  4  2     000 
4   66666       666666666666666666666666666666  4   22    00 
66   666666666b66666666666666666666666666666  44   2?   00 
6   H6   6666666b6666666666666666666b666666666   444  ??  00 
66  BBB    66666666666666666bb666                 44   2?   0 
H  BB   666666666666666666666    444444444    44   ?22   0 
BB  BH   66       66666666666     44         444  44  2?2 
(   BBBB       666666666666     44          444444   2  2????????? 
666    BBM6   666666666b66    444   72?    4444444  ?2  2?   4 

6  BB  BB  6  66666666b6   44   22?  22  4444444   ?   ??  4    444 
666    B   66     6666666  44  7222    2?  222  2?   444 

6666666  666  4  222        2222222227222      2?2?2 

666666  44  2  00      22?     ???? 

44444444444444  22  222    444     72? 

4   2222222       22      0       ?2   4444  44444   ??  2  4  6  4  2   222  22??  ?  44  444 

4   ?         222  22^2     000     22  4             44  2?  44  6  4  ?00         ?2?  44  444   5. 

66                               444444   2     22222222222       222?!       444      6       44  ?   4444   2             ?  44  5 

5                    ,...444444444    722     22      2222     22?    4444       6666666    44     222222222   00     2   4  5 

•  5  6666  4,44,44444,44  ?  22  4  44  22?  2?  44  666666666666666  44  ??  0  ?  2  0000  2?  4  5, 
5  44444444444,4  7?  0  2  44  2222?  44  666666666  6666  44  2  2  4  22  22222  44  5 
5   44                ?2??.???222222          O0      22              444      666666666   ft   666      4  2?  ?2  4444          44  5 

♦6  4   22222      ??7?               000000     222  44444444      66666666666   888   666     44  ?  2   4          444  666  5- 

S  444          4444   7??22??2         0O   2?  44444  6666666   BAB   6666666   4    44  44  5 

5    44444444,     ,4  222                   0  2  4  ,44  444  6666   B  66  666   66  444    b666666666  44   4444     5 

♦5          44   66  66  ,,4,44    22       2  4   6  44  222   44  666   P.  66  66  ft   6    666           66  4   2?   44  6 

444444444444       44                27  44  6   44  44  66       666  ft   66666666666   66   44   2  22   44    6 

3  44444444  4444444,4  44,  4444444444  666  66b66666666  4  72  2  4444  6 
♦222222?   444  22???'.  4,       66    44,444  666     BBBB            66666    44  ?  0  2?  444444 

1       ?22  ?            44444444    444  6b6666666666      444    666    ftriB            666666   4   ?     00   77    444 

1       2222?     4  ??????2222????2  22???   44    b666666       666    444  44   666   BBB             666666  4  2  000000    ??    44 

•1  7   4  44   ??????    00000  7   4    66              66   44   2   44  66   8HBB            66666664  ?  0000000     ?   44 

1  72?     ???»??  22?  2222       444   66             666   44   2   4   66            6666666     444  ?2    0000     22  44 

1  2???      7   4     ?????.  222          4444   66             66   44   22  44  666      666666         44     ???2    00O     ??    44 

*1    22  OOOno  2  4444444   72  2   4444444   66             66   44   22  4  6         666      444444   22  222?  2?  0000    22  44 

44   222  0000O0   7  44     4,4  4,   66606666666   44    2?22  4  6666   666    4444444444   2222722?   00000    2  44 

76  4  22  O000000"0     22222   44  6666  6  4  2   44              4   2222?  2  44   666  666   444    44    44   22?       00000O  2   4 

•  AB6  22  0000000000000000  22  4  6  BB6  42  2?  444444444444  2  7  44  66666  44  2  4  6  44  22  0000  00000  22  4 
76  4  OOOOOOGOOOOOOOO       7  44  666  6   2  2  44              444    22  ?2   4444        4444    4   6   4  22       00      0000  2   4 

4  22  OO0O0  00000OII0  272  44  6  4  2  2  4  66666  66666666  44  2  ?22?  44444444444  6  44  ?22  0  2  4 
•1         00000    U00OUO00              2       4  4   22  4  6     6666          6  4    2722222722   44444      66  4   22  ??  44 

???  222722  22?  2?   7   6   66     6666666    44    ???2?    444    44444   666  4   2???2  0  2  44 

22222222  77?    44   7  00?  46  66  4  44   6   44444     4  ??2   4      6  44   44  6  66  44     2??2     000  2  44 

•??  2222  2?  44,  4  2  ?  0000000  46  64?  44  44  74664?  4  666666  44  2  4  68  66  44444  ?2  0O0O  2  6 
???  2?2  «'•  66  4  2  ??2  222  2  4  66  4  444444  44  6  6  4  ?  4  44  22  4  6  666  4  27722222  2  446 
12?   44444  4444,4  4444,  66  0  42  ??         2??222  44  6   44   4,44444444    6  44  2   44   444     22  4  6    666  4  ?       ?2?2?      4446 

•  12?  4,444  66  44  »  6  ft  4  7?  444444444444444  444444444444  44  ?  4444  ??2?2?2  4  6  64  20  22  444444444 
3  44  44444,4444  66  6  44  44444444444444444444  44  444  4444444  66  6  42  ??2?22  4 
44           44      444     6     4444444444  4444444444    ??2?    ?   4       44444   ????  44        666         6  4   7722  2222772274 

•  666666  4  7?  ??  4  66  44  444  4444444  22  00000  2  4  444444  ?  0?  46  6  44  ?2  1 
h                                 *4    ,,           666  44  ?2    44       4444  2                      2  4  666666666   44  7  7  4   666     666666666   44    22?22?2???21 

5  6*6666    6666  1  ?  22222222?       44  4   ?       2  4  666  64  ?  2  4  66  4  2  44  4444 

•S  666         666  6666   4   1100          24642     22  46    688BBBBB   4  2?  4  6  4   0  2  44         444444444444  44444 

1  666666f.6»  666666  666  4  2  ???222?.?2    44  6   44    2?  ?  4  6            6  4  ?  44  6  4  ?2?   44444   44444444444444   4444 

5  66  66666666666  44              44  666bb66  4  ?  2  4  666b6666666   4    44    44     44444444444444444444         4 

•S  66    66  ,444444  44444   666    BBB   20   4  666         4444444444444444444444444444444444             5 

F  66  666  444,44444444    66600   6  4  2  44  6 

*  666  4444444444  666666     444 
♦S  .'.4  44444444,             66666  44 

5  4  4444444        66        666 

*  4444444  4444444      6666b00 
♦S                            (.4444444        44444444       666     666 

5  444444, 4'., 4       444         6666666 


fi       66666 

4444 

*>.      66666666 

44444 

*>       66666     666 

44444 

S       66       66*6 

4444 

5      66        666 

44 

H      666          66 

44 

5    6666         666 

44 

S   b6666           666 

*   666            6b 

44 

S      66666666666   l*U 

4           6      44  ??????? 

4444444        44   ? 
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FREQUENCY    TAn 

1  F 

InTKR,, 

M 

UNITS 

PERCENTAGE 

CUMUIATIV£ 

NO 

UA1A 

(i 

HELD 

*    LOW 

(1 

0.0000 

0.0000 

.0252 

TO 

1 

1677 

246 

2.0959 

2.0959 

1 .  1877 

TO 

2 

3503 

91  J 

7.7618 

9.B577 

2.3503 

TO 

3 

,5128 

1094 

9.3210 

19.1787 

3.5128 

TO 

4 

,6754 

1348 

11 ,4«50 

30.6637 

4.6754 

TO 

5 

,8379 

2515 

21.4260 

52.0917 

5.8379 

TO 

7 

,0004 

369? 

31  .4561 

83.5478 
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TO 

e 

,1630 

1326 

11.2976 

94.8454 

8.1630 

TO 

9 

,1255 

514 

4.3793 

99.2247 

9.3255 

TO 

10 

,4881 

87 

.7412 

99.9659 

10.4881 

TO 

11 

.651)6 

3 

.0256 

99.9915 

11 .6506 

OH 

OVER 

1 

.0065 

100.0000 

T<)T"I 

1  1737 

100.00U0 

100.0000 

SEARCH  ACTIVITY  -  CiF'-itWATTON  OF  ALTERNATIVES. 


HFPF.ATED  USE  OF  It     MODIFIED  VERSION  OF  THF  HKITISH 
ROAD  RESEARCH  LABORATORY  MINIMUM  PATH  ALGORITHM 
GFNfcRATES  AS  MANY  AS  SFVEN  ALTERNATIVES. 


THE  GENERATE^  ALTERNATIVES 


III                            II 

I  CHOICE  I  PATH  TOTAL  I  RELATIVE  PATH  I  NUMHER  OF  I 
I  NUMHEP  I  I  V>LUE  I  LINKS  I 
ITT                                          11 

ITT  II 
I  1  1  134. 9?S  1  1.000  I  48  I 
III                                             II 

ITT  11 
I  2  I  156. 448  I  1.151  I  44  I 
ITT                                             II 

ITT  II 
I  3  I  162.700  I  1.197  I  48  I 
III                                             II 

ITT  II 
I  4  I  171. /65  I  1.264  I  44  I 
ITT                                             11 

ITT  II 
I  5  1  180.050  1  1.325  1  44  I 
ITT                                             11 
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HAP    OF    GENERA1EU    ALTFOr.ATIVtb. 

«,  choices  joe  sho.n. 
E«CM  IS  LOCatEU  el    a  STRING  OF  nuhrehs  «m 

COHHfSPOND  TO  TMF  CHOICE  NIIHf'EH. 

THE  ALTERNATIVES  arf  OtFINFU  ON  A  NETWORK 
HAVING  ?b    COLUMNS  .*Y  33  ROWS. 
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I)  EXAMPLE      OF  PUNCHED    OUTPUT 

nu"ber  1  of  2  factors»wfight=  1.1)0  earthwork  cist  factor 

number  2  of  2  factors. w£i«ht=  1.00  pavement  construction  cost  factor 
to^htnEo  utility  surface.                                    25  33  2 

5.745  5.736  5.48b  5.610  5.799  6.207  5.792  5.905 

6.307  6,185  6.194  6.248  S.997  6.195  7,681  6.334 

5.675  6.240  5.697  2.210  7.524  1.B20  1.026  1.437 
1  .476 

5.950  5.809  5.722  5.659  6.287  6.201  5.667  5.864 
6.227  6.135  6.095  6.186  5.922  H.053  8.486  6.418 

9.951  6.254  5.31b  2.532  1.924  1.446  .688  1.199 
.997 

5.968  6.212  6.062  5.669  6.443  5,431  6.413  5.620 

6.021  5.478  5.b95  6.470  6.557  9.090  6.767  6.86B 

6.526  5.599  1.162  3.471  8.641  2.962  5.319  1.719 
2.H05 

5.686  5.911  8.103  6.119  8.180  7.803  5.687  6.000 

9.956  10.033  6,242  6.118  6.691  6.698  6.684  6.289 

5.245  4.797  2.825  3.225  7.984  5.437  5.8?3  6.793 
7.128 

8.19*  8.5)6  8.6H5  8.284  6.027  7.745  6.104  7.046 

10.102  6.954  5.486  6.237  6.408  6.366  6.186  5.319 

1.329  2.123  3.186  2.906  1.975  1.616  6.313  6.113 
6.021 

6.230  6.310  6.062  5.553  5.807  5.818  7.832  8.296 

6.290  6.367  6.572  6.439  6.809  6.230  5.7S1  5.268 

1.188  2.628  2.388  7.08b  7.428  6.168  6.483  6.598 
5.823 

6.399  6.032  5.631  5.555  5.701  5.650  5.396  5.438 

5.6bb  5.476  6.014  6.352  5.894  6.448  6.201  6,126 

8.652  10.222  2.074  3.533  6.886  1.4?0  7.701  ,985 
5.595 

6.727  6.657  6.511  6.279  6.014  5.534  5.056  1.806 

5,351  6.625  6.009  6.062  5.598  6.088  5.797  9.750 

3.356  2.640  3.805  6.030  2.518  5.154  2.152  5.580 
5.539 

6.921  7.060  8.970  6.637  5.762  5.465  4.853  5.694 

5.855  6.322  6.487  6.393  5.861  5.019  3.450  .747 

2. [58  6.666  4.964  .697  6.187  6,179  5.830  6.315 
5.56? 

6.608  6.8^?  9.054  8.670  6.435  5.6P?  2.349  2.076 

6.09b  6.508  6.493  5.886  .688  6.46?  2.887  4.068 

8.016  6.055  2.024  5.346  7.250  6.416  7.007  6.266 
6.012 

6.069  8.064  8.6?9  8.349  6.074  5.750  5.4?6  5.753 

6.356  6.673  6,450  5.607  8,333  4.545  8.298  2.172 

?,889  1.753  .940  5.275  6.430  6.657  6.516  6.383 
5.966 

5.681  5.660  6.315  5.561  .575  2.339  6.247  6.599 

6.771  6.994  7,873  7.727  7.687  7.671  7.2?2  1.879 

.991  1.513  2.025  3.188  6.356  6.485  6.318  6.116 
S.H49 

5.47?  5.484  5.859  ,234  5.043  6.164  9.814  7.665 

7.345  7.663  7.111  7.312  7.367  7.146  7.00?  4.820 

1.067  1.388  6.292  5.941  6.547  6.345  5.861  5.927 
5  ,666 

6.292  .025  6.184  6.231  6.156  6.974  11.102  7.947 

8.491  7.862  7.394  6.870  4.387  3,948  6.166  2.676 

1.221  2.156  5.709  5.537  6.182  5.366  5.292  5.416 
5.497 

?.222  5.648  9.225  6.353  6.777  6.770  6.974  11.357 

6.940  7.361  7.316  4.477  1.879  5.774  5.0^0  1.392 

6.417  5.304  4.550  2.137  1.889  7.7?5  1.337  6.045 
5.778 

h.Mo  D.981  7.618  9.1V9  6.905  6.595  7.442  6.797 

6.440  6.582  6.832  1.534  1.694  1.048  1.884  4.257 

2.606  1.439  5,306  ?.8?4  7.948  1.746  6.440  6.256 
5.50? 

7.339  6.213  7.284  6.907  6.803  6.411  6,0?1  6.*?7 

1.657  1.527  3. 5^3  2.29?  ,830  1.797  6.054  6.848 
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7.224 
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5.256 
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2. 235 
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6.449 
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6.256 

6.869 

5.  '16 

5.290 
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5.470 

5.21  1 

1  .630 
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5.4?4 

.487 

1.958 
6.955 
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9.096 

3.659 
6.737 
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7.119 

.891 
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.794 
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2.151 
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1.593 

6,976 
3.050 
4.559 
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9.154 
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4.137 
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